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Abstract : We compute the one-loop coefficients for an alternative Symanzik improved pure 
gauge SU(iV) lattice action (N = 2 and N — 3). For the standard Symanzik improved 
action we confirm previous results by Liischer and Weisz. 


1 Introduction 

When studying pure lattice gauge theory it is important to realize that the lattice spacing 
a is a rather elusive quantity. The action, when formulated in terms of the natural link 
variables U^x), does not contain any a dependence. The reason is that the continuum 
action is scale invariant. The introduction of a lattice spacing by parametrizing 

U^(x) = exp(aA /i (x)) (1) 

and insisting that A M is the (dimensionful) continuum vector potential is only natural in a 
classical context where one can limit oneself to smooth fields A^. In a quantum mechanical 
context, however, defined by the path integral (with coupling constant go) 

Z = JvUex V (-S Lat [U]/g 2 0 ), (2) 

typical field configurations are wildly fluctuating at the scale of one lattice spacing, and 
eq. (|T|) is merely a convenient way to make contact with other regularizations of the con¬ 
tinuum theory. The introduction of A M is not at all necessary for the extraction of physical 
quantities, for example by means of a Monte Carlo simulation. 

Nevertheless, the notion of a lattice spacing is sensible in the quantum mechanical 
context too. The reason is that eq. (|2|) dynamically produces a length scale, namely the 
correlation length £. More generally, it produces masses of particles, m*. Of course due 
to scale invariance only dimensionless ratios appear, fhi = arrii (■fhi being the mass in 
lattice units). This means that the continuum limit is defined by m* —> 0. Moreover, the 
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proximity to the continuum can be assessed by considering ratios m^/rrii as a function of 
a. For small a (i.e. small fhj) one thus expands 


Ulj Tflj ( 2 .') 9 ( 2 ) 4 

—— — — Rest 27 5 Rcstjfl — Caa a -I - c ^ 7 - cl -)-•••, 
m l m l J J 


( 3 ) 


where the coefficients are independent^ of a. Note that in the presence of external 
physical mass scales, for example as set by a finite volume L 3 , one can of course also 
consider ratios with 2ti/L. 

For the conventional Wilson action |lj all coefficients c-J generically are non-zero. 
Improvement means reducing the lattice artefacts (‘Rest’) for all values of i and j that are 
associated with physical masses. By Wilson’s renormalization group argument 0 it should 
be possible to fold a lattice action that produces Rest = 0. In a certain approximation, 
the implementation of this approach has been put on a practical level by Hasenfratz and 
Niedermayer et al. [3!]. In this paper we follow the less ambitious perturbative approach 
due to Symanzik Q. 

Symanzik improvement amounts to expanding ‘Rest’ with respect to a, and systemat¬ 
ically removing the a 1 2 , a 4 , • • • terms. Usually only the removal of the a 2 term is pursued, 
and we will do so in this paper. From now on we mean by Symanzik improvement this 
limited reduction only. For pure gauge theories, Symanzik improvement is believed to be 
attainable ||| 6 | by adding to the standard Wilson action a finite number of extra Wilson 
loops, for example^: 



Due to asymptotic freedom 0, the continuum limit a —> 0 is achieved by approaching 
the critical point, go —■► 0. Therefore it is sensible to calculate Cj(gg) perturbatively as a 
function of g 2 } . To this end one expands 


c,(g 2 0 ) 


£ c, W (9„ 2 ) 


m=0 


( 5 ) 


Symanzik improvement to n loops amounts to fixing c® • • • c 4 - n ' ) such that physical quanti¬ 
ties show only 0(a 4 , a 2 g^ n+1 ^) deviations from their continuum values. At tree level (n = 0 ) 
the analysis is rather easy, and many alternative Symanzik-improved actions have been in¬ 
troduced. However, for a long time the availability of a one-loop Symanzik-improved action 
was limited to the Luscher-Weisz action [| 8 |. In ref. [|9| we have presented the one-loop co¬ 
efficients for a second action, namely the square action that was introduced in ref. | 10 ] . 
The purpose of the present paper is to give the details of our computation. 

In present-day Monte Carlo simulations g$ is typically of order 1, and this value seems to 
be too large for Symanzik improvement to be effective. For this reason a phenomenological 


1 For convenience we neglect logarithmic corrections for the moment. 

2 Taking the liberty to deviate from the notation in ref. [[| we assign C 4 to the 2x2 Wilson loop. 
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extension of Symanzik improvement, called tadpole improvement, was introduced [|TTf . 
In this approach a mean-field parameter |]T2[ u 0 ~ {U^x)) is included in a tree-level or 
one-loop Symanzik-improved action by dividing each link variable by uq. This has been 
argued Bi to capture most of the higher-order corrections in go, thus approximating an 
all-loop Symanzik-improved action. In practice, Monte Carlo data for tadpole-improved 
actions often, but not always, show dramatic improvement. As a side remark, let us 
mention that the goal of an exact all-loop Symanzik action might be attainable to sufficient 
precision by using Monte Carlo techniques. In ref. |L4| this has been studied for the 
fermionic sector of QCD, where the leading lattice artefacts to be cancelled are of the 
order a instead of a 2 . 

We believe that our introduction of an alternative one-loop Symanzik-improved action 
is useful for several reasons. It allows one to make the following analyses of improvement: 

• Estimate remaining lattice artefacts by comparing the effectiveness of two different 
Symanzik-improved actions. 

• Study the universality and consistency of tadpole improvement. 

(At tree-level such analyses were performed in, e.g., ref. jnj). Moreover, a one-loop calcu¬ 
lation has some spin-offs: 

• The Lambda parameter ratio A improve d/Awiison, which is needed for doing simulations 
with both actions at the same physical scale. 

• The tadpole parameter u 0 to one-loop order. 

The outline of this paper is as follows. In section we discuss on-shell improvement [6j, 
which implies that only two improvement conditions need to be satisfied in order to cancel 
the 0(a 2 ) errors in all spectral quantities. In section 0 we set up lattice perturbation 
theory. All propagators and vertices needed are listed in the appendices 0-0. In sectiong 
we calculate the static quark potential to one-loop order. This allows us to extract the 
Lambda and tadpole parameters, and also one improvement condition. In section |6] we 
discuss perturbation theory in a twisted finite volume, which subsequently is used to obtain 
two improvement conditions and, again, the value of the Lambda parameter. As an extra 
check on the universality of on-shell improvement we compute one additional spectral 
quantity. We conclude with section [7j, where the main results are listed in a compact way. 
For the sake of readability we start in section with discussing the physical setting of 
the twisted finite volume computation. We will do so from a somewhat more conventional 
perspective than was given by Liischer and Weisz in refs. 
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Our calculation follows closely the well-documented work of Liischer, Weisz and Wohlert 
Therefore we only go into full detail at points where we take a different 
approach, and at some points where the earlier treatments can be complemented. Special 
attention will be paid to the many checks we have performed to convince ourselves of the 
correctness of the final results. 


2 Twisted boundary conditions and electric flux 

In this section we review twisted boundary conditions in SU(iV) pure gauge theories that 
were introduced by ’t Hooft |L7|, and apply the analysis to the setting of section |6|. For 
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convenience we restrict ourselves to the continuum formulation. However, twist is also 
well defined on the lattice [I 8 |, and the generalization of the results below to the lattice 
formulation [|B| used in section [| is straightforward. 

Twisted boundary conditions are defined as follows (no summation convention): 


AJx + L„0) = A-(x), 


( 6 ) 


where z> is the unit vector in the positive x v direction. L v denotes the size of the (rectan¬ 
gular) space-time in the v direction, and A^ v is a gauge transform of A^\ 


A“(x) = fl(x)(^(x) + d^Q-^x), O(x) G SU(iV). 


A structural requirement is: 

^ii(^ + T^r / )hl i /(x) -|- 12^(x )z 


(z = e 2wi/N 


), 


( 7 ) 


( 8 ) 


for each /i, v plane. This is a consistency relation, necessary because there are two ways 
to relate x to x + L^fi + L u v. n^ = —n u/1 are integers, only defined modulo N (i.e. 
n lw e 7L n ). This guarantees that Z fW = exp(27r in^/N) is a center element of SU(iV). In 
pure gauge theories it is not required that = 1 , because A^{x) is invariant under gauge 
transformations with a center element. If = 0 (rriodAQ, it is possible to render the twist 
matrices O,: ^-independent by a gauge transformation. Below we will assume that this has 
been done. 


How pure gauge theory is quantized in the presence of twist is described in refs. ||T 7 | . |TT) 
The outcome (neglecting instanton effects 


is that physical Hilbert space associated 
with a twist splits up in N 3 sectors, which are eigenspaces of gauge transformations 
compatible with the boundary conditions. The Hamiltonian can be diagonalized separately 
in each of these sectors (which are usually labeled by a vector e € In this context an 

important operator is the Poyakov line 


Pj(x^) = Tr 


Pexp 


dxjAj(x.) 


a 


( 9 ) 


(x (1 ) = (; x 2 ,x 3 ), and cyclic). Under a gauge transformation compatible with the boundary 
conditions it picks up a phase factor (e Zjjv), and from this it can be shown that it maps 
the e th electric sector into the (e + j) th sector. 

Twisted boundary conditions thus allow for N 6 quantum sectors: N 3 due to (physi¬ 
cally) periodic boundary conditions in space, labeled by , and within each of these: N 3 
sectors labeled by e. These sectors have a beautiful interpretation, due to’t Hooft |T7|: e 
corresponds to electric flux, m to magnetic flux (with = — Yhk £ ijk rn 'k)- For shortness we 
will not repeat’t Hooft’s arguments (see also ref. [l||) why this is a natural identification. 

We would now like to apply the above results to the specific geometry used in section j6|. 
To this end we restrict ourselves to a magnetic flux m = (0,0,—1). We also choose 
Li = L 2 = L and send L 3 —> oo, thus obtaining a ‘twisted tube’ (this terminology is 
due to Liischer and Weisz). It can be proven that the vacuum valley of the twisted tube 
consists of only one point, A^ = 0 (up to a gauge transformation), which is why it admits 


a straightforward perturbative analysis, contrary to the purely periodic case m = 0 [21 . 
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What we want to do in the remainder of this section, is to classify the lowest-lying 
excitations above the quantum vacuum |0). To this end we use a path integral approach. 
We can get a handle on the spectrum by considering the gauge invariant correlation function 
of two Polyakov lines at a distance t in (Euclidean) time, 

(Pl(x 2 ,x 3 ;t)P l (x2,x 3 -,0)). (10) 


Since |0) is a state with zero electric flux, this correlation function is associated with the 
creation of a (gauge invariant, i.e. physical) state at time 0, with electric flux e = (1, 0, 0), 
and its annihilation at time t. Therefore for large values of t it will be dominated by the 
exponential behavior exp(— Mt), where M is the energy of the lowest-lying Hamiltonian 
eigenstate that has non-zero overlap with A(^ 2 ,^ 3 )| 0 ). 

In order to compute m perturbatively, the allowed fields in the path integral must be 
parametrized. To this end one defines the x-independent matrices 

T[n] = H- n2 H , 2 l W5 (ni+n2)(?ll+n2 - 1) , n G TL\. (11) 


The choice of the phase factor is a matter of convenience. The important property of T[n] 
is that it satisfies = z Uv (v = 1, 2). From the second paper in ref. |18| we know 

that 

^Tr (r[n]TV]) =<W (mod A0- ( 12 ) 


This allows us to Fourier transform A^x) in the following way: 


a ( r ) — 

^ ’ (2tt) 2 L 2 N 


E / dk 0 dk 3 e ikx A,(k)T k . 

k u k 2 Jr 


(13) 


Here kx = Y^=o k^x^ and J2k u k 2 runs over 


2tt 

k v = mn v , m = ——, n v G 7L , (u — 1, 2). 

1 V 1j 


(14) 


Also we adopted the notation of ref. [15 


Tfc = T[(ni,n 2 ,0)], {k v — mn v , v = 1,2). 


(15) 


In contrast to the purely periodic case, A^(k) G (D carries no SU(1V) index. Instead, the 
momenta in the twisted directions are quantized in units of m rather than 2^r/L = Nm. 
The total number of degrees of freedom is preserved, because Tr A^x) = 0 implies 

Afj,(k) = 0 for 7ii i2 = 0 (mod N). (16) 


Using this formalism one quickly derives (with t = Xq) 

Pi(x 2 ,x 3 ;t) = f dkodhe^A^k) +0{A\), 

\ ZT[ ) L ng^ JE “ 

k = (fco, 0, (2ir/L)n + m, k 3 ). 
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(17) 

(18) 














Note that k 2 is now quantized in units of Nm, but carries an additional momentum m. 
Physically, this is the Poynting vector e x m due to the electric flux e = (1,0,0) created 
by the Polyakov line ||22]. 


For the moment we neglect the 0(A\) corrections in eq. (|T7]). We then obtain 

(Pi (x 2 , x 3 , t)Pi(x 2 , x 3 ; 0)) = / r2 dk 0 dk 3 e lkot D^\k), (19) 


nGZS 


where D^ l (k) is the dressed propagator, defined through the 2-point function: 

(A;(k)A u (k f )) = ^((2n) 2 L 2 N)8 n ' ijni 5 n ' 2in2 8(k' 0 - k 0 )5(k' 3 - h) Xk D™(k). 


( 20 ) 


The factor = (1 - $ nu o( mo dN)8 n2 ,o(modN)) is implied by eq. ([16D : modes with zero mo¬ 
mentum (mod (Nm)) in the twisted directions are absent on the twisted tube. In eq. (f20P, 
/q, = mriy and k! v = mn! v (is = 1,2). In eq. (|H?D, the momentum k is restricted to eq. (|T8|). 
Now it is well known that 


lim f dk 0 e ikot D[f(k) ~ e~ E{ ^\ 
Jr 


( 21 ) 


where k = (LE(k);k) is the pole of D^(k). To lowest order in perturbation theory this 
gives the usual formula E(k) = |k| = \Jk'l + k% + k 3 . Here it is essential that the Polyakov 
line, due to its gauge invariance, couples only to physical (i.e. transversal) polarizations. 
Indeed, from eq. (|T7[) we see that P\(x 2 ,x 3 ) couples to the polarization = <5^1, which is 
transversal due to eq. (|T8D . 

Inserting this result in eq. (|H|) we see that the leading exponential decay comes from 
k = (0, m, 0). Thus the lightest particle with electric flux (1, 0, 0) has a mass m (+C?((7g)). 
This is what Liischer and Weisz called an A meson. From eqs. (|14|), (|l^) we see that it is 
the lightest particle on the twisted tube. 

There are more particles with the same mass. To lowest order in perturbation theory 
these are associated with the other poles of the bare propagator for |k| = m, i.e. at k = 
(m, 0, 0); (0, —m, 0); (— m, 0, 0). Of course they are created by Polyakov lines with electric 
fluxes e = (0,—1,0); (—1,0,0); (0,1,0). It can be proven that the total of four particles 
now found have equal masses to all orders in go. This is simply a consequence of symmetry 
considerations. Namely, the Hamiltonian has two CP symmetries, one corresponding to 
(1) Xi —> —Xi, and the other corresponding to (2) X\ <->• x 2 . We will not prove this here 
(for the proof see ref. |L5!). Let us only mention that (for N > 3) the charge conjugation 
C (defined by A^x) —> HL4*(x)IF -1 , where W G SU(1V) can be freely chosen) is not 
a symmetry, because it violates the boundary conditions (it changes m into — m). The 
behavior of the electric fluxes under the two CP transformations (referred to as CP 1 - 1 ’ 2 ^) is 
given in table [L|. 

One may wonder about the O(Af) corrections in eq. (^). In principle they bring about 
O(gg) corrections to the above analysis, because they give rise to 3-point (or in general 
n-point, n > 3) functions in the expansion of ([[(]). The leading Feynman diagram of this 
type is given by (in coordinate space) 

ti 
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(ei,e 2 ) 

CPW( ei ,e 2 ) 

CPW( ei ,e 2 ) 

(1.0) 

(1,0) 

(0,-1) 

(0,1) 

(0,-1) 

(-1,0) 

(-1,0) 

(-1,0) 

(0,1) 

(0,-1) 

(0,1) 

(1,0) 


Table 1 : Properties of the electric fluxes created by the Polyakov lines Pi, P2, P\ and P 2, 
under two CP transformations; o; 2 ) = (—aq,a; 2 ) and ~P^ 2 \xi,X 2 ) = ( X 2 ,X \). 


For t —> 00 , we must have t\/t —> 0. The reason is that from the above analysis we 
know that each propagator of length t' suppresses the correlation function by a factor 
of approximately exp(— mt'). Tims only diagrams in which a single (dressed) propagator 
runs from the one Polyakov line to the other contribute to the leading behavior for t —> 00 : 
exp(—(m+(P((yfg))t). Diagrams like the one above do not contribute to the O(g^) corrections 
to the mass m. Instead, their contributions should be associated with the renormalization 
of the composite operator Pi(x 2 ,x 3 ). 

There is a further set of four particles which at tree-level are degenerate with the four 
particles found above. These also carry electric flux |e| = 1, but they have a different 
polarization, = 5 Mi3 , at momenta k = (0, m, 0); (m, 0, 0); (0, —m, 0); (— m, 0, 0). We will 
refer to this set as A - mesons, because under appropriate CP transformations they have 
eigenvalues opposite to their counterparts found above. Those we will call A + mesons. 
Due to the CP (1,2 ) symmetries, all A - mesons are exactly degenerate. However, there is 
no symmetry relating the A - mesons to the A + mesons, so beyond tree-level they must 
be expected to have different masses. The A - mesons can be created by gauge-invariant 
operators, but these are of a more complicated nature than the Polyakov lines creating the 
A + mesons. 

No more particles with masses m + O^g^) exist, because pure gauge theory admits only 
two physical polarizations. From the bare propagator one sees that the next particles^ 
encountered in the spectrum have twisted momenta \ki\ = |fc 2 — m, and thus their energies 
equal y/2m + (P(gg). These were called B mesons by Liischer and Weisz. In our language 
they carry electric flux \ei\ = |e 2 | = 1. Still higher up are particles with masses 2 m + O^gfy 
or more. For <70 7^ 0 (or rather for the renormalized coupling gn^L) 7 ^ 0) many of these 
particles should be expected to be unstable, because the kinematics allow them to decay 
into A or B mesons (only the particles with tree-level masses 2m might be safe, depending 
on the O(gl) corrections). 


3 On-shell improvement 

It is important to specify precisely what quantities are to be improved. I 11 Symanzik's 
work [Q the aim was to improve correlation functions. However, correlation functions 
depend on the elementary field operators. As a result it is difficult to find a formulation 
in terms of which this kind of improvement can be implemented. For the 0(3) a- model, 

3 It is natural to consider particles with momenta (0, m, fc 3 ), /c 3 7 ^ 0, as A mesons in motion. 
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Symanzik succeeded by performing a well-chosen field redefinition, but for gauge theories 
the problem becomes more intractable and to our knowledge no one has solved it up to 
date. 

One can avoid this complication by restricting oneself to ‘on-shell improvement’ as 
introduced by Luscher and Weisz |J. Instead of improving correlation functions one now 
pursues only the improvement of spectral quantities. Unlike correlation functions these are 
necessarily invariant under field redefinitions. Examples of spectral quantities are masses 
of stable particles and the static quark potential. 

In order to determine a minimal set of improvement conditions it is sufficient [6j to 
perform a small-a expansion of the classical action. The most general action we consider 
is given by eq. (f3|), the expansion of which reads (see e.g. refs. [H§HD : 


SdQ^o)}) = —y(c 0 + 8ci + 8c 2 + 16c 3 + 16c 4 )(po) 51 Tr (f^x)) 

+ yx(c 0 + 20ci - 4c 2 + 4c 3 + 64c 4 )(^q) Tr ( D ti F ^{ x )) 

+j(c 2 + 3c 3 )(g 2 ) Y, Tt(d,F, x (x)D u F uX (x)) 

X,[l,V,\ 

+jC 2 ( So 2 ) y Tr (D^fxjj + 0(0*). (23) 

x,fi,iy, X 


In this equation, F^ v = 


8 u A fl + [Ap, A u \ and D M — 8^ + adA M , where A M is related 


to U^ by a path ordered exponentiation 
appears only in the combinations 


see 


eq. (|247|) ). A crucial observation is that c 4 (#q) 


co(^o) = c o{gl ) - 16c 4 (c/o), Ci(^q) = ci(^o) + 4 c A (gl). (24) 

This implies that with the substitutions Co(go) —*► co(5 , o)> ci(fi'o) ~^ Ci (go), the entire 
analysis of ref. || can be copied literally. The only subtlety is that positivity of the lattice 
action must be rechecked for c 4 (#q) ^ 0. Postponing this proof to appendix [E|, we take 
over the following conclusions from ref. |J: 

• The normalization condition 


c o(do) + 8ci (^o) + & c 2(9o) + 16c 3 (<7o) — 1 (25) 

can be imposed. This ensures that eq. ( |23| ) has a natural continuum limit and, more 
importantly, determines Co(go) as a function of the other coefficients. Since in the 
path integral the action is weighted by a factor of l/<?o, eq. ( |25| ) is always attainable 
by a redefinition of the bare coupling constant go. 

• After imposing the normalization condition, on-shcll improvement is expected to be 
satisfied for all spectral quantities once the coefficients c\(gl), c 2 (go) and c 3 (go) are 
tuned correctly. Furthermore, c 3 (gg) can be chosen 0 because of the freedom of field 
redefinitions. Hence in fact only c 4 (gg) and c 2 (go) are to be determined, i.e. we need 
only two improvement conditions. 



• The tree-level conditions for on-shell improvement are 

1 

Cl ~ _ 12 ’ 

c 2 = 0 . (26) 

Moreover, the fact that co(gij), c\ (c/q) and c 4 (go) enter the analysis only through the 
combinations Co(. 9 o) an d Ci(c/g), hnplics that one of them, e.g. c 4 (gg), is a free function. 
Having noted this, one should yet keep in mind that the one-loop coefficients^] c' Q , c[, and c ' 2 
do depend on the tree-level choice of c 4 . In particular, the Liischer-Weisz and square actions 
will have different one-loop coefficients. Nonetheless c' 4 is a free parameter for both actions. 

In the remainder of this paper we limit ourselves to actions satisfying eq. (p5|) and 
C 3 (g'g ) = 0, though some intermediate results are valid more generally. In particular we 
focus on the Wilson action [l| (abbreviation W; c 0 = 1, Ci = c 2 = c 3 = c 4 = 0) as well 
as the Liischer-Weisz || (LW; c 0 = §, c 4 = — c 2 = c 3 = c 4 = 0) and square JTIJ (sq; 
c o = T; c i = — i, c 2 = C 3 = 0, c 4 = 134 ) Symanzik actions. For the last two actions we 
determine the coefficients c[ and c' 2 in sections |] and |^. 


4 Lattice perturbation theory 


This will be a very short description of lattice perturbation theory. For a thorough setup 
of the general formalism see refs. [|U||24 . 


The elementary field in pure lattice gauge theories is the link variable U tx {x). In a 
perturbative approach it is convenient to use the parametrization ©■ Since U^(x) G 
SU(IV), A^x) lies in the Lie algebra of SU(IV): 




(27) 


where H“(x) G 1R. We limit ourselves to the fundamental representation, with conventions 


Tr (r a T b ^j 


= -hS, 


2 ^ ab1 


T-.T- =f dx T c . 


(28) 


Note that eq. (|1|) contains an ordinary exponential, not a path ordered exponential. There¬ 
fore eq. ([h|) receives corrections. In appendix |Q this provides us with a non-trivial check 
on the vertices. 

Using eq. ([[]) it is straightforward to expand the action ([|) to a given order in g 0 . The 
only complication is the many terms involved, but this can be handled by a symbolic 
computer language such as FORM [^5|] or Mathematica |26| . 


Another issue is gauge fixing. We have adopted the gauge choice of ref. [|I 6 

ntage of 

0. 


i.e. 


covariant gauge fixing. The advantage of this choice over the Coulomb-like choice in ref. |15 
is that the gauge condition, 

, covariant g.f. 


(29) 


4 We use the notation d = c[°\ c[ = c, 


(i) 
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is independent of the lattice action chosen, thus eliminating a potential source of mis¬ 
takes when changing from the LW action to the square action. Here d L is the left lattice 
derivative, d^f(x) = ( f(x ) — f(x — afi))/a. We implement eq. by adding to the action 

S g{ = ~ a4 J2 Tr > ( 30 ) 

a X \ fj, ) 


where a is the well-known gauge parameter that should drop out of physical quantities. 

Of course the gauge fixing procedure induces ghosts in the usual way. Hence a ghost 
term S g h os t must be added to the action. One more term to be added is the Haar measure 
term ^measure associated with the parametrization (JT|), (|2j]). The total action thus reads 


9 ■-’'total 

9o 


—2 (S({Ci(go)}) + S g {) + ^measure 
% 


+ S, 


ghost • 


(31) 


For the calculations in the following sections we need the expansion of S tot& \/ up to third 
order in g 0 . It is given in the appendices 0-0 Using this expansion, and the Euclidean 
path integral 

Z = f exp , (32) 

perturbation theory is set up like in the continuum, with the exception that coordinates are 
now discrete or, equivalently, momenta run over one Brillouin zone only. Correspondingly, 
momentum-conserving delta functions 5{k\ + • • • + k n ) are 2n/a periodic. In the remaining 
sections we use lattice units. The a dependence can be reobtained from a dimensional 
analysis. 


5 Static quark potential 

5.1 Generalities 

This section is devoted to the calculation of the static quark potential V(L) to one-loop 
order. Its definition reads 

V(L) = Jim -i In (i (Re TcU(L x T)>) . (33) 

Here U(L x T) denotes a rectangular Wilson loop with extension T in the time direction 
(0) and L in some spatial direction, say 1. The expectation value (• • •) denotes averaging 
with respect to the path integral in infinite volume. Since V (. L ) captures the exponential 
decay of a gauge invariant correlation function, it can in principle be expressed as an 
(L dependent) function of the spectrum. Therefore its small-a behavior, or in lattice units 
large-L behavior, is suitable for a (partial) determination of the on-shell improvement 
coefficients. 

To avoid confusion, note that we will not pursue the computation of V(L) in the 
confining regime, i.e. at large physical distances L £ (where £ is the correlation length). 
That would not be feasible within the framework of ordinary perturbation theory. However, 
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it is the basic postulate of Symanzik improvement that all physical quantities can be 
improved simultaneously, so that considering short-distance quantities only is sufficient for 
a complete determination of the improvement coefficients. 

The potential V(L) not only depends on L , but also on the coupling g 0 . Roughly fol¬ 
lowing the notation of refs. 00, but limiting ourselves to the fundamental representation 
of the Wilson loop, we expand 

/I \ 00 n 2n 

- In (- (ReTrf/(L x T))J = ^jL- Wn (L,T). (34) 

We restrict ourselves to the tree-level and one-loop contributions (n = 1,2). In terms of 
Feynman diagrams they read 



The endpoints of the gluon lines ( x/x/x/v ) represent source termsQ due to the Wilson 
loop, and hence are to be summed over all lattice sites on the loop. On the other hand, 
vertices ( • ) should be summed over all sites in space-time. The striped graph represents 
the one-loop vacuum polarization: 



The first diagram on the right hand side represents c[ insertions (cf. eq. (|180|) ), the second 
one stands for the measure term insertion, eq. (ED- The other symbols have the usual 
meaning, with ghost propagator ---->— . As in ref. 0] we have split the four-vertex in 
two parts, corresponding to the second and third lines of eq. (|183|) . 

The analytic expressions corresponding to the graphs in eqs. (|35|)-(|37l) can be found in 
ref. [0], both for finite T and for T —> oo, and in terms of arbitrary coefficient^ c*. Below 
we copy those expressions whenever needed. 

At tree level the following result was found p|,[Tb|]: 


lirn T) = 4C/ [ sin 2 (i/qL)A)o(fc)k=o, (38) 

_ T^o o 1 _ J k 

5 These are proportional to go, cf. eqs. (Jl|), (|27|). 

6 In refs. ||[L(|, C 2 <-> C 3 with respect to the convention of refs. ii® and this paper. 
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where 



Cf is the Casimir in the fundamental representation, 



( 39 ) 


(40) 


and D^ u (k) is the gauge field propagator, see eq. (|181| ) and appendix |Bj. Since /i 
and k 0 = 0 , the form Q206| ) of the propagator is applicable: 


D oo( k ) k ~° k 2 - (ci - c 2 )k (4) - c 2 (k 2 ) 2 

= k 2 - (0! - c 2 + ^ 4 ) - c 2 (k 2 ) 2 + 0(k 6 ). (41) 


The symbol k (4 ) and similarly k^ 4 ) are dehned in eq. (|197|) , however now summing over 
spatial indices only. By doing a contour integration one finds [p) (up to an irrelevant 
L-independent term) 


T lto iu,(i,T) = -TL 


1 + 3 


Cl — C2 + J2 

U 2 


+ 0(L~ 4 ) 


(42) 


which obviously is improved for c± — c 2 = — 

Now let us proceed to the one-loop level. In ref. [[Ip] it is shown that 

lim l- w 2 (L,T ) = 48C/ f sin 2 (±/ciL)-u; 2 (k) (D 00 (k) | fco=0 ) 2 , 

T —>oo i ' J k 


(43) 


where the L-independent function w 2 ( k) has the expansion^] 

f „ k^ 4 ) ~ 1 

w 2 (k) = k 2 < ai — /3 0 lnk 2 + (a 2 + b 2 lnk 2 )^- + (a 3 + 6 3 lnk 2 )k 2 + 0(k 4 ) > . (44) 


The coefficients a„ 


are zeroj if eq. ([26|) 


b n depend on the tree-level values ci and c 2 . In particular, b 2 and b 3 
holds, which from now on we will assume. On the other hand, 


A) 


11 N 
3(47t ) 2 


(45) 


takes precisely the universal value required for renormalization (see eq. (|50|) below). In 
order to express the improvement condition in a clean way we perform the integration in 
eq. (|43| ). It can be shown that S 3 does not contribute to the L dependence of w 2 , and the 
result reads (up to an irrelevant constant) 


1 fiO 

lim -w 2 (L,T) = —L 

T^oo i 7 TL 


Ahr - 


a 2 / a 4 N 

Si + (3qc — 3a 2 — + O ( — 


(46) 


7 We adopt the convention that 0(a n ) may also stand for a n In a terms. Note that in lattice units, a is 
absorbed in k. 

8 This is to be expected from Symanzik’s work on scalar theories 0. For the present situation Weisz 
and Wohlert flu] have shown it explicitly for C 4 = 0. Extension of their result to other tree-level improved 


actions is not difficult, as we will see in subsection 5.2 


12 



















where c is a numerical constant (c = —2q = —1.1544313298 ■ • •). For clarity we reinstated 
the a dependence. Note that for tree-level non-improved lattice actions eq. ([R]) receives 
additional 0(a 2 /L 2 ) corrections. 

The total one-loop static quark potential has a finite regularization-independent limit 
for a —> 0, 


V(L) 


Ci_ 

AttL 


9 2 r 


1 + 9r 


3ci2T7 + 0 




(47) 


in terms of the renormalized, L-dependent coupling 


9l = 9l 1 - 


1 % hi ( — ) + A,c - a± 


9 2 o + O(g 4 0 ) 


(48) 


Since V(L ) and gn are a-independent physical objects, we deduce from eq. 
loop improvement condition 

= 0 . 


the one- 


(49) 


The calculation of the c[ dependent coefficient 02 is described in the next subsection. 


Also Lambda parameter ratios can be extracted |]16|.h7| [29|] . As g r must be independent 
of both a and the regularization prescription, it follows from eq. (^) that go depends on 
a and cl, C 2 . This is expressed by the well-known renormalization group flow 


a —>0 


Po9o 


ln(a 2 A 2 ) + (!?(ln — ln(a 2 A 2 ) ) , 


(50) 


where the Lambda parameter A depends on the lattice action: 

= e 5 fe(ai (Li })-ai ({<*})) 


A 


(51) 


In this equation {c*} and {c*} represent any two choices of tree-level coefficients that satisfy 


the normalization condition (|25|) , and A, A* are the corresponding Lambda parameters. 
As a final remark, also Wilson loops for finite values of L and T are objects of interest. 


They allow for a gauge invariant definition of the mean field or tadpole parameter uq [JTTJ]: 


Uo= ( — (ReTrC7(l X 1)) 


(52) 


More generally one might define for L, T e IN 

/ 1 \ 2 (L+T) 

u 0 (L,T) = y— (ReTr f/(L x T))j , (53) 

but the choice L = T = 1 is probably optimal. The reason is that if L or T becomes larger 
(i.e. closer to the correlation length), Uo(L, T ) may start picking up physical signals instead 
of the lattice artefacts that tadpole improvement aims to correct for. Nevertheless we will 
also compute tfi(l, 2) and Wi(2, 2) because for the LW and square actions they are related 
to the vacuum polarization. 
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Liischer-Weisz 

square 

/(1,1) 

0.366262680(2) 

0.3587838551(1) 

h 1.2) 

0.662626785(2) 

0.6542934512(1) 

n (2,2) 

1.098143594(2) 

1.0887235337(1) 


Table 2: Several numerical values of eq. (|55|) for the Liischer-Weisz and square actions. 
Note that H (L, T) = /(T, L). 


5.2 Particulars 


We continue with the details of the one-loop calculation, ft should be kept in mind that 
we will only compare the Wilson, Liischer-Weisz and square actions, and therefore we set 
to zero all tree-level coefficients except Co, c\ and C4, and all one-loop coefficients except 
Cg, cl, c' 2 and c' 4 . 

As an appetizer let us start with finite Wilson loops. From ref. [|l(J we read 


with 




W\(L,T) — C f /(L,T), 


r f sin ^k\L\ 2 /sin ^k 0 T\ 
*k l sin \k\ ) l sin | k 0 I 


Dioio(k)- 


D[iv\ P is the propagator of —i(d a A T — d r A a ), 

= { y.MV. -(/*<-* ")] - [A « p 


(54) 


(55) 


(56) 


and J k is the four-dimensional analogue of eq. (|5D|). Upon insertion of the expressions given 
in appendix [B|, the evaluation of /(L, T) for hnite values of L and T can be handled to a 
high accuracy by standard integration routines (we used NAG pi] for all numerics in this 
subsection). The results for L, T e {1,2} are listed in table ^]. The entries (L,T) = (1,1) 
and (L, T) = (1, 2) for the LW action agree with the results of Weisz and Wohlert. They 
found /(l, 1) = 0.366262 and /(l, 2) = 0.662624. 

The rest of this subsection is devoted to W2(k), defined in eq. (|43|). Eqs. (|36D and (|37|) 
carry over in 


W 2 ( k) = U^(k)+U^(k)+7T 00 (k), 


n 


= 7 t ' + 7T 


meas + ^ghl + 7r gh2 + ^ + ^ + ^W 


(57) 

(58) 


where k = (0,k). The corresponding analytic expressions are copied below from ref. JT6 


Propagators and vertices can be looked up in the appendices. For shortness we use the 
notation D = D 00 and adopt the summation convention. 

U {i \k) = ND~\k) [ (2n) 4 5(k^- k'+ k") \ (± - lD-\k)D(k")) 

Jk',k" I ' ' 

( 59 ) 
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( 60 ) 


t/ (3 >(fc) = iND-'lk) / (27r) 4 <5(fc + k' + k") 

Jk',k" 

Bt Dfi0 (k') Du0 ( k ") Vo , u (k, v, e 


71 


V 

meas 

/ii/ 

7T ghl 
n fiv 


(k) 
C k ) 
(*) 


(2 k' 0 

(M, 




_j_ /VA 

12 1 v U /XI/J 
24 


k„5 


liV\v 


) q'ii\(k)kx, 


7T 


gh2 

/it/ 


(*0 = 

J k 


\N I (27r) 4 6(k + k' + k") 

k',k" 


k,K- (k 1 - k") (k' - k") l 


7r 


£(*o 


■IN 


7r 


7T 


(*) = 


(*) = 


k' 2 k" 2 

f (2n) 4 S(k + k' + k")D\\i (k')D pp , {k") 

Jk',k" 

xV pXp (k, k',k")V u y P '(k, k',k"), 

$N [ D Xp {k') [V Xppu (k', -k', k, -k) - V Xppu (k', k, —k', -k )\, 

J k' 

k/N^) 'k> V,, f K^(k', -k',k,-k)D^x(k'). 


(61) 

(62) 

(63) 

(64) 


(65) 

( 66 ) 

(67) 


In eq. (^TJ) q' denotes q pu , eq. (|194|) , with, following the prescription of appendix |A]. tree- 
level coefficients c* replaced by their one-loop counterparts c'. Note that because eq. (|25|) 
is satisfied at tree level (and 03(5^) = 0), 


I6C4. 


Cq = — 8c^ — 8c' 2 

The tensors Kfy in eq. (|67| ) are dehned in appendix |C]. The unusual factor 

d = ^(6 C, -N) = (iN - Ji) 


( 68 ) 


(69) 


arises from a contraction [JT6[] of eq. ( |186|) . 

It is worthwhile to analyze how the above expressions lead to eq. (fb|). In the first place, 
they are completely finite for non-zero k because UV divergences are regularized by the 
lattice (integrations run from —n to 7r) and potential IR divergences are regularized by the 
external momentum k. For k = 0 some of the integrals are logarithmically IR divergent, 
and this carries over into the In k 2 corrections! to the polynomial behavior in eq. (f44|) as 
a function of k (we come back to this point below). As mentioned in appendix [A], even 
(odd) vertices are parity even (odd). This also holds for the two-vertex, i.e. the propagator. 
From this the invariance of eqs. (|59|)-(|67|) under k —» — k can be verified explicitly. Hence 
■u)2 (k) contains no odd powers in k. Also cubic invariance can be checked, and therefore 
up to order k 4 only the combinations k 2 and k^ 4 ) appear. Finally, there is no constant 
term because 7T pi/ (0) = 0 holds [0 due to gauge invariance while the expressions for tA 3,4 ) 
clearly are of order k 2 due to an overall factor of D~ l (k). 

For use below note that 7T Aljy (0) = 0 allows us to subtract from each term on the right 
hand side of eq. (|58|) its value for k — 0. 


9 In non-lattice units this gives rise to the familiar In a UV divergences. 
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We now embark on the computation of the various contributions to u> 2 (k). Note that 
all of them are proportional to N, except for ti' 00 and We therefore follow Weisz and 
Wohlert and define 


*( k ) = Ko( k )\k 0 =o-,x(k) = nZ(k)\ ko=0 -y( k) = i(n) 2 (k) - x(k) - z(k)). (70) 

x(k) and z(k) are easily computable, 

x(k) = d {Jc 0 /(1,1) + 8ci B (1, 2) + 16 c 4 H (2, 2)j k 2 

- Cl /(l,2) + 4c 4 /(2,2)]k( 4 )}, (71) 

*(k) = (ci - c 2 )k (4) + c 2 (k 2 )^. (72) 

The elaborate part of the computation is posed by y(k). We can save ourselves a fair 
amount of work by pursuing only the following two objectives: 


1. The evaluation of A sq ( uare )/A W ( ilson ). Due to eq. 
where di is defined in eq. 


this requires computing cq q — of , 


2. The evaluation of the improvement condition (^), i.e. of the coefficient a 2 defined 
in eq. (Ill), for the square action. It is sufficient to compute a 2 q — a 2 w , because the 


value of a 2 w was computed by Weisz and Wohlert 
and Weisz fOl . 


, and later verified by Liischer 


We thus need only consider y*(k)—y**(k), where V and denote different lattice actions. 
This implies that we can drop vr^, eas , ttoq 1 and 7r qq 2 . Much more importantly, it gives us the 
opportunity to circumvent the logarithmic terms in eq. (|Tij ) due to the following argument. 

The most straightforward way to perform the k —» 0 expansion of the integrals would 
be to expand the integrands. Unfortunately this method is too naive for terms that contain 
both propagators D(k') and D(k"), because due to momentum conservation the expansion 
of D{k") gives for each factor k 2 a factor ( k ,2 + C>(/d 4 )) _1 , thus possibly introducing IR 
divergences. As an example consider n^l(k) (eq. (|65|)). For k = 0 its integrand behaves as 
1 jk' 2 near k' = 0, so the integral is convergent. However, if we expand the integrand with 
respect to k 2 , we introduce at 0{k 2 ) a logarithmic IR divergence. This would be fatal were 
it not that due to its IR nature the divergence is completely determined by the small-U 
behavior of the integrand. In particular it is the same for all lattice actions satisfying 
eq. (|25|) . This implies that subtracting the integrands of two such actions produces an 
extra factor of ( k' 2 + 0(k' 4 )) +1 , and the 0{k 2 ) term is finite. 

For other contributions to ?/(k) one can give similar arguments. Therefore the naive 
method does yield the correct k 2 -coefficients in the expansion of ?/*(k) — y**(k). Further¬ 
more, if V and ! **’ both denote tree-level improved actions, then the integrands are the 
same up to srifeleading order, and the naive method also produces the correct /^-coefficients. 

As a side remark we mention that if one is interested in y( k) itself, instead of merely its 
Ci dependence, one has to tackle the expansion with respect to k in a more sophisticated 
way. For example, one might subtract a suitable integrand that has the same IR behavior, 
but the integral of which is computable by continuum methods. This leads to the Ink 2 
factors appearing in eq. (|44l) . From the above discussion we then conclude that the coeffi¬ 


cient {3q is universal within the class of all lattice actions satisfying the condition (|25p, and 
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the coefficients fe 2 ,3 are universal within the subclass of tree-level improved actions. Using 
Weisz and Wohlert’s work [ TTj ] on the Liischer-Weisz action it then follows that the latter 
values are b 2 = 0 and 63 = 0. 

Let us denote by Na * the contribution of Ny(k) to a*. From eq. ( fUD and the above 
discussion we conclude that for the computation of of 1 — it is sufficient to consider 
( 7r oo( fc ) - ^00(0))/^. froo ( k ) - *oo(°))/W, and 


m 3 \k) = il? [ ^^D0(k')Drt(k')V o ^(O,k',-k'), 
Jk ' (2fcf' 

- *jA\- 


k ' (2A& 

1 

6 (2 V Q ) dk' 0/ 


cos \k' Q d 


D(k’). 


(73) 

(74) 


Similarly, a 2 q — a^ w only receives contributions from — 7:^(0))/N and (^00 (^) — 

-Kqq(0))/N (C/( 3 > 4 ) (/j) only produce k A terms of the type (k 2 ) 2 ). 

The integrands of IT 3 ’ 4 ) are already expanded with respect to k. The expansion of 7 Tqq (k) 
can be performed rather easily, because the single internal propagator is independent of 
the external momentum k. We used Mathematica m to obtain 

1 (ffiTo( fc ) -7roo(°)) 


N 


k 2 

12 




+Cl 


Co + 8c x ) [(7 - fk' 2 ) D(k') + k r 0 k[D 01 (k') 

(84 - 27k' 0 2 - 60k' 2 + 10£q 4 + f k'fk'f) D(k') 


k( 4 ) 


12 


+2 (18-7 k' 2 ) k' 0 k[D 01 (k') 

+c 4 [2 (-60 + 23k' 2 ) (3 - k' 2 ) k' 2 D(k') 

+8 (—12 + 5ky) (—3 + k 2 ) k'^k^Doi^k') | 

|ci (25 — 6 /cq 2 — 10 k' 2 + 2£)o 2 &q 2 + D{k! ) 

+ (4-A:; 2 ) k' 0 k[D 01 (k') 

+c 4 [(l2 - 3k' 0 2 - 120k' 2 + 38A)o 2 A)i 2 + 38Aq 4 - f k' 0 2 k[ 4 ) D(k') 
+2 (12 - 5k' 2 ) (4 - k' 2 ) hfoDo^k ’)]} . 


(75) 


For c 4 = 0 this agrees with ref. [16 


The last ingredient needed is the expansion of (n^(k) — ttqq(0))/N. In principle the 
analytic expansion of the integrand is straightforward, but to our taste too tedious due 
to the fc-dependence of the propagator(s). We therefore decided to evaluate A7 t^’*(A;) = 


( n oo’ Sq (k) ~ <r(0) — 7Too ,l {k) + 7Too’*(0))/iV (i = W, LW) numerically for a number of 
values of k. As explained above the result should be of the form 

At T^’\k) = Aa^k 2 + Aa^ 3,i k< 4) + Aa[ 3li (k 2 ) 2 + 0(k 6 ). (76) 

(For i — W there are also logarithmic terms at (T(/c 4 )). We determined A a( 3 ’ w by choosing 
k = (e, 0, 0) and fitting to e 2 . For the determination of Aa^ 3,LW we subtracted results for 
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k = (e, 0, 0) and k = (e/ a/2, ejs/2, 0), thus eliminating the k 2 and (k 2 ) 2 terms. The results 
were fitted to \e A . These methods introduce numerical instabilities of the order of 8/e 2 
and 8/e A respectively, where 8 is the computer precision (<5 rs 10 -14 ). One would like to 
choose £ not too big, but on the other hand subleading terms (which are of relative order 
£ 2 ) must be kept sufficiently small. By chosing £ in the ranges 0.0001 to 0.001 and 0.01 to 
0.1 respectively, we were able to extract A a/ 3 and An)]" 3 with roughly the same precision 
as we expect we could have attained by an analytic expansion and subsequent numerical 
integration. 

Our final result in the y(k) sector is 

a 8 / -af = -0.031810197(2), 

a s 2 q -4 w = 0.000087063(4). (77) 

As a check of our analysis and programs we also computed 

4 W - af = -0.031361443(2), (78) 


which agrees with the value found by Weisz and Wohlert jl6|, -0.03136145(1). When we 
include the results for x(k) and z(k), using k 2 = k 2 — ^k (4 ^ + 0(k 6 ) and table |2] (together 
with /(l, 1)Wilson = §), we obtain for the coefficients of u>2(k) 


~sq 

a 2 


~ a r = 


-LW 


"—(0.086580342 (3) A — 0.082828348(1)^) for i = LW 


(o.085608020(2) A— 0.08069673318(6)^) for i = sq, 


-21 


sq 


^2JLW 


0.002158351(4)A - 0.0033681212(1) 


A 


(79) 


(80) 


It follows from eq. (|5lD that 

^lw 

A w 
A. 


A' 


sq 

W 


4.1308935(3) 

5.2921038(3) 

4.0919901(2) 

5.2089503(2) 


for A = 2 
for A = 3, 

for A = 2 
for A = 3. 


(81) 


The values for the LW action agree with those in ref. [IT]: 4.13089(1) for A = 2 and 


5.29210(1) for A = 3. Moreover, the Lambda ratios for the square action completely agree 
with the results in ref. ||3T|. These were obtained in a totally different way, namely by using 
a background field approach. 

From eqs. (|49|), (|80|) we read off a one-loop improvement condition. We can substitute 


the result of Luscher and Weisz [15] for the LW action, 


lA ^Jlw 


—0.01100879(1) 
—0.02080086(2) 


for A = 2 
for A = 3, 


(82) 


which is consistent with the value extracted by Weisz and Wohlert [16|] from the static 
quark potential, but more accurate. In this way we obtain 


A - ^2Jsq 


—0.00837615(2) 
—0.01544851 (3) 


for A = 2 
for A = 3. 


(83) 
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We conclude with the observation that the last diagram in eq. (|37|) gives quantitatively 
far the most important contribution to these values: 



n w 

[ci - c' 2 

]lw 

- a~ 2 sq 

\tt w 


0.78 for N = 2 

0.82 for N = 3, 

0.69 for N = 2 

0.75 for N = 3. 


(84) 


The other diagrams give much smaller contributions, so that the large ratios are not due 
to coincidental cancellations. This fits well in the tadpole/mean held picture because the 
PU-vertex is not present in the continuum theory, and hence is at least partly responsible 
for the deviation from 1 of the tadpole parameter uq. 


6 Spectroscopy in a twisted finite volume 


6.1 Introduction and formalism 


In this subsection we give a rather short overview of the formalism of pure lattice gauge 
theory in a partly twisted space-time. For details we refer to refs. H, [L5|. and also to 


section [2]. As mentioned in section we limit ourselves to infinite volume in the Xq, £3 
directions, and twist in the £ 1 , £ 2 directions (with twist quantum — 1 ): 


A^x + Lu) = n^A^x)^ 1 , (v = 1,2), 
where G SU (IV) are A M and x independent matrices satisfying 

0 , 0 2 = z = e 2 * i,N . 

Eq. ( |T3|) generalizes to 

A„(x) = Jo$e“ ( " +i “' i) T( fc )ri., 

k 

where the integration symbol ‘$’ stands for 


r n dk 0 r^dkz 
k L 2 N ^ ’ 




and T k G SU (N) satisfies 


QUtW 1 = e ik - L T k , (f = 1,2). 


(85) 

( 86 ) 

(87) 

( 88 ) 

(89) 


The momentum components k± = (fci, fc 2 ) in eq. (| 88 D are discretized as in eq. ( 0 ). but on 
the lattice J2k ± runs over, say, n v — 0,1, • • •, NL — 1. The solution to eq. (|89|) is given by 
eqs. ( 0 ), ©• 

I 11 terms of 


(k, k') 
(k, k') 


( 0 if n v — 0 (modV), (v = 1, 2) 

1 1 otherwise, 

n-in\ + n 2 n 2 + ( n i + n 2 )(n[ + n 2 ), 

nin' 2 — 7i\7i 2 , 


(90) 

(91) 
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the matrices satisfy JTj| 


IV = 

r, 

k if Xk'-k = 0 

T fc = 

1 

X 

0 

Tr T fc 

= 

0 

t-b 

X 

?r 

rt — 

1 fc — 



r fe r k 1 


: r fc+fc /Z5< fc . fc '>-3( fc . fe ') 


( 92 ) 


From these properties, eq. (^7]), and the fact that A^(x) is in the Lie algebra of SU(1V) it 
follows that 


Mky = -ztevA^-k), ( 93 ) 

A^k) = 0 if Xk = 0. (94) 


For later use it is convenient to define a delta function associated with eq. (|88|) : 


S(k) = NL 2 5 k0 2ir8\k 0 ) 2irS\k 3 ), (95) 

where 5 1 is the ordinary 1-dimensional (27T-periodic) delta function, and 5 kk , is the NL- 
periodic 2-dimensional Kronecker delta, 



1 if n v — n' v = 0 (mod(ATL)), (^ = 1,2) 
0 otherwise. 


(96) 


Eq. ([7|) is an extremely nice property of the twisted tube, at least in perturbation 
theory. It implies that A M (0) is not a degree of freedom of the twisted gauge held, and 
therefore is not to be summed over in loop ‘integrals’. In this way the twisted tube escapes 
any IR problems that lure in infinite or finite periodic volumes. 

A strongly related effect can be seen very clearly at tree level. As follows from the 
Feynman rules listed in appendix 0, the bare propagator equals 

{A lt (k)A v {k , )) go=0 = 5{k + k') (-§*“*(*-*)) XkD^(k). (97) 


(We slightly disagree with ref. Hl5 
Feynman propagator D^ v (k) and thus reads 


Its mass-shell equation is determined by the ordinary 


k 0 = ±iE„( k), S„(k) = jk? ± + kl (1 + 0(k 2 )). 


(98) 


However, since k± is discretized and k± = 0 is not permissible, a mass gap emerges, of 
width m (in the continuum limitf 10 }) . As long as L is chosen so small (in physical units) 
that asymptotic freedom ensures g 2 R (L) -C 1, the tree-level value of the mass gap receives 
only small quantum mechanical corrections, which are computable in perturbation theory. 

In refs p|,p~5|, Liischer and Weisz present a Kaluza-Klcin picture of the twisted tube 


by viewing it as a two-dimensional theory; the compact dimensions are considered internal 
space. In this interpretation rii and n 2 are quantum numbers. The particles defined by 
n_ i_ = (1,0) and n± = (1,1) (or cubic transformations thereof) are called ‘A’ and ‘B’ 


10 Exact lattice formulas can be found in appendix 
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mesons. Their masses are = m and m-g = \[2 m, to lowest order in the lattice spacing 
and coupling constant. The mesons possess a ‘spin’ quantum number taking values ±, 
corresponding to the two physical polarizations of the underlying four-dimensional theory 
of massless gauge bosons. For an understanding of these particles in terms of the electric 
fluxes introduced by’t Hooft, see section f|. 

The twisted tube brings about many well-defined spectral quantities that can be used 
to extract the Symanzik coefficients. The simplest quantities are the A and B masses. In 
subsection |6.2| we compute to one-loop order the mass of the A meson with positive spin. 
We will see that this gives the same improvement condition as was found from the static 
quark potential. The reason is that in both cases the coefficients c[ enter via a propagator 
insertion. We nevertheless decided to undertake the calculation in order to become familiar 
with perturbation theory on the twisted tube, and also to check eq. (|83l). 

The Symanzik coefficients c\ and c' 2 can be obtained separately by also computing 
the one-loop correction to the three-gluon vertex. For suitable external momenta this 
is a spectral quantity, because it is related to the clastic scattering amplitude of two A 
mesons || (which in turn is related to the energy eigenvalue of a two-particle state). The 
calculation is described in subsection |6T . 


6.2 Mass of the A + meson 


The energy spectrum of (improved) lattice field theories is best defined through the trans¬ 
fer matrix J52J. The spectrum dictates the exponential decay of the two-point function 
(A p (x)A p (y)) as a function of the time separation xq — r/o- Equivalently, the spectrum can 
be read off from the pole structure of D 1 ^ 1 (k), the definition of which in the twisted tube 
geometry we already encountered in eq. ( f20| ) (cf. eq. (0)): 


(A^k)A w (k')) = S(k + k‘) (-$*-**■*>) joDjfft). (99) 

We concentrate on the A meson with positive spin. From the discussion in section [| 
we know that to all orders in g o it can be represented by the polarization e p = 8^\ at 
a momentum k = (fc 0 ,k); k = (0,m, k 3 ). Let us denote the corresponding eigenvalue 
of D™\k) by d {ull (k) (= D^\k)). The energy E( k) and the associated wave function 
renormalization Z (k) belonging to this A + meson are defined through 


dfuii(fc) 


Zfr) 

kl + E 2 (k) 


+ (regular in k 0 ), 


( 100 ) 


where the expansion is valid for /c 0 close to A'(k). We will assume that A(k) is physical in 
the sense that A(k) ~ a 0 instead of ~ a -1 , in non-lattice units. This can be verified from 
our results below. It implies [|32| E( k) e 1R. For future use we note 


Z~\ k) 


J__d_ 

2 k n dkc 


v/ -1 

“full 


(k) 


k 0 =±iE(k) 


( 101 ) 


In terms of the vacuum polarization defined through 

D™\k) = D^ik) + glD, p {k)E pa {k)D av {k) + C% 0 2 ), (102) 
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we find 

4diO) = d (o)( k ) + 9o d lr) (k) + 0{g £), d^(k) = -tt u (k). 

For the one-loop energy and wave-function renormalization, defined through 

E( k) = E^+glEM + O^), 

k) = Zo(k)+^Z 1 (k) + 0(^ 0 4 ), 

the following formulas are valid: 

1 


(103) 


£i(k) = 


2 Er 


— Z 0 (k) d^(fc) 

(k) uv 7 dA > ko =iE 0 { k )’ 


Zl(k) = -z 2 (k)< 


1 d 


1 d 


- 2 Bo ( k )£ l( k, ^U ( - 0 ;w 


(104) 

(105) 

(106) 


1 k 0 =iE 0 (k) 

For the mass calculation we put /c 3 = 0, since m A + = E(k = (0, m, 0)). At tree level, 


eq. (|206|) and hence eq. ( |211|) are applicable, so that [|| 

= m — (c'i — c 2 + ^)m 2 + 0(m 4 ). 


(107) 


Note that this is improved for the choice (|26|) of coefficients, as it should be. Below we 
assume eq. (^6| ) is satisfied. 

At one-loop level we write 


m A + = m { A+ + 9 o m A+ + O(go), 
so that from the above analysis it follows that 


mjjj = -Z 0 ( k) 


TTll (k) 


2 m 


( 0 ) 

A+ 


(108) 


(109) 




Note that this quantity only depends on m and (implicitly) N, or equivalently L and N. 
From eq. ( |211|) we find 

^(k)l t . ( „, m , 0) = l + O(m 4 ), (110) 

For the extraction of the one-loop coefficients we can neglect the 0(m 4 ) correction term. 

The vacuum polarization n fJ , u (k) has of course the same expansion as in section [5|, see 
eqs. (^7[) and fl58P , but the explicit forms ('Gp-(^) change a little on the twisted tube (as 
follows from appendix 0): 


^u(k) = - (k\5» v ~ k^xu) q'f, x (k)kx, 

7'2 

^\k) = ±5^$g 2 _(k,k')jE, 

k' k - 

i q ^ , „ _ kl v ~ (k^-k"), 

= -J$ S(k + k J ' + k")£(k,k')-^—'- 


k',k" 


2 U"2 


k' 2 k 


ini) 

( 112 ) 

(113) 

(114) 
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( 115 ) 


= i $ S(k + k' + k")g 2 _(k,k')D xx ,(k')D pp .(k") 

k',k" 

k', k"), 

rft( k ) = -*$9-(k,k')D Xp (k') 

b' 


= 

x W^xpiAv^k i k , /c, fc) , /c, k , A))], 

1 (k\8, lv - k^Sx^kx °i $ (Xk' + k ')) 

i k' 

(116) 


y.Kfl(k',-k',k,-k)D^ x {k'). 

(117) 

In these equations 

g-(k-k') = 2isin fe')) , 

(118) 


which equals 0 if Xk 1 = 0 or Xk+w = 0. A reasonable check of the above formulas is that 
they equal their infinite-volume counterparts in leading order for m —> 0 (i.e. L/a —> oo 
in non-lattice units), which is what one would expect physically. The technical reason 
is that in this limit the factor g 2 _(k,k / ) is oscillating infinitely faster than the rest of the 
integrands, so that it can be replaced by its average value: —2. Due to the smoothness of 
the resulting integrand, $*,/ can be replaced by N f k , and %k'Xk' by (TV - TV -1 ) f k ,. 

In ref. [|T5| it is shown that near the continuum limit m^j is of the form 


m 


(i) 

A+ 


m 


~ (i) ~ (i) 

—— — (c[ — Cn) m 2 + 0{m A ); —— = a 0 + aim 2 + 0(m 4 ), 

m m 


(119) 


where we separated the contribution of 7(containing all dependence on c(). The coef¬ 
ficient^] a* are determined by the other contributions. It is clear that the improvement 
condition reads 

c[ — c 2 = a\. (120) 

In eq. (|119|) there are no quadratically (1/m 2 ) or logarithmically (m° lnm) divergent terms 
because m^ J { is independent of go, while multiplicative renormalization of go alone is suffi¬ 
cient to cancel all divergences in the continuum limit. Also the m 2 lnm term is absent due 
to tree-level improvement, cf. the discussion below eq. (]T^). Note however that individual 
diagrams can give 1/m 2 , lnm or m 2 lnm contributions to m^|/m. For example, 7 r meas 
contributes N/(2Am 2 ). In the Coulomb gauge, odd powers of m can even appear [|15|, but 


not so in the covariant gauge. The reason is that the Coulomb propagator contains 1/k 2 
divergences for k —> 0, while the covariant propagator only has 1/A: 2 poles. In any case, 
?h^ + /m being of the form ( |119|) is a very good global check against computational errors. 

Following ref. [T3], we decided not to perform the small-m expansion of eqs. ( |113| )-( |117| ) 
analytically. Instead, we computed the sum of Feynman diagrams numerically for a number 
of values of L, and fitted the results to the expected form ( |119| ). This is much alike to what 
we did in the previous section to compute the coefficients appearing in eq. m- However, 
a difference of practical importance is that the CPU time needed for the evaluation of 
m^+ increases for decreasing m (see below), while for A7 Tqo(A:) it is independent of the fit 

11 In spite of the same notation, these are different from the coefficients in subsection |5.2|. 
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variable k. Another difference is that + depends on N. Due to the twist it is difficult to 


unravel this N dependence analytically, and like in ref. [|T5| we did separate computations 
for N = 2 and N = 3. 

We did not use entirely the same approach as in ref. [15 . 

are: 

1. We used covariant instead of Coulomb gauge fixing. 


The most important differences 


2. Liischer and Weisz fully automized the generation of vertex subprograms. That 
is, they wrote meta-programs that take a Wilson loop as input, and give vertex 
subprograms as output. We wrote the subprograms by hand, using the vertices that 
we had already used before for the computation of the static quark potential. As 
mentioned in section f|. the generation of the vertices themselves was automized by 
means of Mathematica |26|] . As a precaution against programming errors, we always 


3. 


performed a number of numerical tests of our subprograms against the corresponding 
Mathematica representations. 

A further technical difference is that Liischer and Weisz used PL/I [Q for the gen¬ 
eration of vertices. 

Liischer and Weisz followed a rather unconventional way to extract the coefficients clq 
and (i\ . For example, to find a 0 they constructed from the data for, say, L — 6, 8, ■ ■ • 30 
a new data series for L — 8,10, ■ ■ • 28 that is improved in the sense that the m 2 term 
in eq. (|119|) is cancelled, leaving only 0(m 4 ) deviations. This procedure was then 
iterated. We preferred doing a least-squares fit, which for our data appeared to give 
somewhat more stable results. 

Once the vertex (and propagator) subprograms are ready, the programming of eqs. 
(|113| )- (|117|) is easy. The main point of interest is the integration routine, since it illustrates 
once more the merit of the twisted tube. One should appreciate that the integration routine 
is required to be extraordinary accurate, because many digits are lost when extracting the 
coefficients, in particular a\. We required a relative accuracy of 10 -13 or better. Whether 
or not an integration routine is capable of attaining such a high accuracy in a reasonable 
amount of time depends on the the number of integration variables and the smoothness 
of the integrand. Now for the twisted tube the ‘integration’ symbol $, eq. (|88|) , involves 
only two integration variables, ko and k 3 . The components k\ and ki are to be summed 
over, costing a rather cheapQ factor of (IV 2 — 1)L 2 . More importantly, the integrand is 
periodic (over one Brillouin zone) and analytic on the domain of integration of /c 0 and 
k 3 (due to the mass gap, poles are shifted into the complex plane). Such a situation is 
ideal for constructing efficient integration routines, especially if one performs a change of 
variables [15] that shifts the pole further away from the domain of integration. As a result, 
we typically needed only 50 2 points to approximate f kl k2 within the required accuracy 
(Liischer and Weisz quote 32 rather than 50, apparently due to the different pole structure 
in the Coulomb gauge, or a more efficient change of variables). 

Concerning the analyticity of the integrand in ko and k 3 , we stress that the implementa¬ 
tion of momentum conservation for diagrams with two internal propagators is to be chosen 

12 Nevertheless this factor is, together with stability considerations, the reason why very small values of 
m cannot be reached. 
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Liischer-Weisz 

square 

L 

N=2 

N=3 

N=2 

N=3 

4 

-0.021583156919197 

-0.040661797276715 

-0.020085367798214 

-0.038849497567066 

6 

-0.019599952046337 

-0.040131915718096 

-0.018908108548525 

-0.039452155788701 

8 

-0.018461969616743 

-0.039274516206705 

-0.018063236142647 

-0.038897070150935 

10 

-0.017889223811184 

-0.038810590019211 

-0.017632105742522 

-0.038571321147693 

12 

-0.017569982136115 

-0.038547675035293 

-0.017390811874792 

-0.038382470120619 

14 

-0.017375041457146 

-0.038385977285131 

-0.017243147194999 

-0.038265025623147 

16 

-0.017247595811504 

-0.038279865552201 

-0.017146489209537 

-0.038187471819923 

18 

-0.017159817570292 

-0.038206617036467 

-0.017079864357098 

-0.038133728057528 

20 

-0.017096835258186 

-0.038153983851014 

-0.017032035100709 

-0.038095009442688 

22 

-0.017050132787137 

-0.038114916592240 

-0.016996555876247 

-0.038066217556170 

24 

-0.017014554076415 

-0.038085133354648 

-0.016969519951262 

-0.038044238228378 

26 

-0.016986831369867 

-0.038061914071922 




Table 3: Raw data for m£j/m. For each entry the absolute accuracy is estimated to be 

(fo 2 ■ 10- 14 . 


wisely. The reason is as follows. Due to tree-level on-shcllness, the external component k 0 
of the momentum k = (&o,0, m, 0) is imaginary, and approximately equal to im. Hence if 
one took k' 0 , k' 3 6 R as the integration variables, the k" propagator would be singular for 
{k' 0 cs k' 3 0, k[ — k' 2 — — m} since in that case k" ~ 0,0) and k 2 ~ 0. It is 

better to shift k' 0 = k 0 — \k 0 (hence k({ = — k 0 — |/c 0 ), where ko (together with k' 3 ) is the 
IR-valued integration variable. 

Our results for individual lattice sizes are summarized in table [|. Note that the can¬ 
cellation of quadratic divergences decreases the accuracy by a factor of 1/m 2 (i.e. l/(am) 2 
in non-lattice units). As mentioned above, we used least-squares fits to determine the co¬ 
efficients a 0 and cq. We checked that the data is consistent with the absence of odd terms 
in m, and also with the absence of terms In m and m 2 In m. We then fitted the data to 
the form ao + aim 2 + 02m 4 + &2'ar 4 hi m + 03m 6 + b 3 m 6 lnm. Our fits strongly suggest that 
62 = 0, but we have been unable to find a rigorous analytic proof for this. 

We found the most accurate results by using a minimal set of points for the fit, 
L = 16 ■ • • 26 for the Liischer-Weisz action, and L = 14 • • • 24 for the square action. The 
disadvantage of this approach is that the error estimate of the coefficients is necessarily ad 
hoc. As error estimate of a* (i = 0,1) we typically used a* — d t , where a; was obtained 
by dropping either 03 or 63 as a fit parameter. The justification for this procedure is that 
that the coefficients 03 and b 3 themselves can barely be determined from our data. As an a 
posteriori justification, the results below show the error estimate to be realistic whenever 
a comparison can be made to other results. 

Using this method we obtain 



0.0168265791(7) for N = 2 
0.0379274963(15) for N = 3, 

(121) 

0.01100890(15) for A^ = 2 
0.0208015(5) for IV = 3, 

(122) 
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a 


sq 

0 


a 


sq 

1 


-0.0168265790(15) for IV = 2 
-0.037927497(2) for N = 3, 

-0.0083763(3) for N = 2 
-0.0154489(7) for N = 3. 


(123) 

(124) 


Note that the continuum coefficient a o is independent of the action chosen, as it should 
beQ. Moreover, using eq. (|120|) we see that the values for a^ w are consistent with the 


somewhat more accurate valued obtained by Luscher and Weisz |15| (these are copied in 
eq. (|82P). The values for a] q agree with our static quark results, eq. (|S3|). 

We would like to report that for the LW action we have also done a computation of 
(which couples to 7733(^0; 0, m, 0)), for N — 3. As expected, we found a continuum 
coefficient different from ao(A" 1 "), namely ao(A _ ) = 0.0000679470(6). A structural check 
on Symanzik improvement is that can be improved simultaneously with (and 
the static quark potential), as we found a^ w (A“) = —0.0208011(3). For N = 2 we did a 
run for the Wilson action, and found a 0 (A“) = 0.00006029(14). Our results for ^(A 1 * 1 ) 
agree with those obtained from dimensional regularization by van Baal 
the analytic spin-splitting formulas 


35|, which imply 


a 0 (A + ) - a 0 (A ) 


1 

67r 2 
3 

87T 2 


for N 
for N 


2 

3. 


(125) 


6.3 Effective coupling constant 

In ref. || Luscher and Weisz define an effective coupling constant A through 

\J Z (k)Z(p)Z(q) ^2 r 3 ( ^ , l-,p,2\q,j) = i\f{k,p,q). (126) 

3 = 1 

A number of new symbols show up. The external lines ( k , 1) and ( p , 2) correspond to 
on-shell A particles with positive spin, while (q, e) corresponds to an on-shell B particle 
with positive spin: 


k = (iE{ k),k); 

k = (0, m, ir ), 

(127) 

P = (~iE{ p),p); 

p = (m, 0, ir), 

(128) 

"cr 

o' 

q = (—771, —m, —2 ir), 

(129) 

e — (0; 1, —l, 0). 


(130) 


Note that these particles and polarizations are completely physical, as they can be created 
by Polyakov lines, see section From now on k, p and q always have the above special 
meaning. The energy E and wave function renormalization Z were defined in the previous 

13 The discrepancy between our result for ao(N = 3) and the value quoted by Luscher and Weisz is due 
to a misprint Q in ref. JE|]. 

14 Liischer and Weisz gained two digits by doing a mass calculation in three compact dimensions. Un¬ 
fortunately this setting is unsuitable for the computation of the other improvement relation by use of 
scattering theory. 
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subsection^. The value of r is defined by E'(q) = 0 (we choose r > 0). Note that at tree 
level, eq. (j206|) applies to all external propagators, with = <5 Mj1 , 5 Mj2 or e M . In particular 


(131) 


A(k) — r — \\plm to leading order in m and go- Furthermore the equalities 

£(p) = £(k), Z(p) = Z(k) 


hold to all orders in perturbation theory, because k, p are related by symmetries of the dis¬ 
cretized twisted tube Jl5|. r 3 (fc, 1; p, 2; q,j ) is the three-point function (A 1 (k)A 2 (p)Aj(q)), 
but dropping the trivial factor S(k + p + q) and amputating the external lines in the usual 
way. Finally, the color factor f(k,p,q) is defined in appendix |Dj. The coupling A is a 
suitable parameter for on-shell improvement, because its square is proportional to the 
residue of the pole in the scattering amplitude for (A + , A + ) —> (A + , A + ), appearing due to 
B + exchange Q. 

We expand 


a = Jo {a> 0 > + Jo 2 a< 1 ) + o( 9 *)}. 


(132) 


The computation of the tree-level value A® simply amounts to substituting the momenta 
k, p and q in the relevant expressions in appendices BO and B (only hs 2) is not listed, 
but it can be looked up in ref. [[L6|). The result reads 


Z 0 { k) = 1 - (ci - c 2 + is) m 2 + 0(m 4 ), 

Z 0 { q) = 1 + (ci - c 2 ) m 2 + 0(m 4 ), 

J2 e i r T e (k, l;p,2\q,j) = -2f(k,p,q)g 0 Y j e j V 12j {k,p,q) 

3 3 

= if(k,p, q) x {-8g 0 m) [l- (4ci- 3c 2 + ^)m 2 + 0(m 4 ) 

It immediately follows that [|] 

A( 0) = -8m {l — \m 2 [9 (ci - c 2 + i) + 2c 2 ] + C>(m 4 )|, 

so that c 2 = 0 is indeed the second on-shell improvement condition at tree-level (eq. 

From now on we assume c 2 = 0 and c± — c 2 = — 

At one-loop level the following formula can be derived from eq. (|126| ) and the expressions 
given in the previous subsection (which can easily be extended to the B + particle): 


(133) 

(134) 


(135) 

(136) 

D)- 


AW 


m 


(l - ^m 2 ) 
-2 (l - & 


4 d 


p(i) 
m k 0 dko 
d 2 


TTll (k) 


ko=iE( k) 


dql 


\Y, e i e 3^iM) 


1,3 


+ CA(m 4 ), 


(137) 


90=0 


where TW is defined through 


X':,r.,!fc. l;l>, 2 : q, )) = l%fi. k p,q) {r 1 " 1 + 9 ,)r : l1 + Ofa’i j. 


(138) 


3 =1 


15 Momenta implicitly label the particle type, and hence no additional label for the Z factors is necessary. 
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In the above equations we may use tree-level expressions for the external momenta, i.e. we 
are free to redefine 


k = (i£? 0 (k),k); 

k = (0, m, ir 0 ), 

(139) 

p = (-*^o(p),p); 

p = (m, 0, ir 0 ), 

(140) 

q = (0,q); 

q = (— m, —m, —2 ir 0 ), 

(141) 


where Tq > 0 is the solution to E'o(q) = 0. This redefinition brings about only corrections 
to eq. (|132|) . One should however be aware of the fact that in the derivation of eq. (|136| ) 
we used eqs. (|139| )- (|T4T|) rather than eqs. (|1 27| )—( |129|) . A priori this can cause corrections 
to A*- 1 ). However, we have checked that all such corrections are at least of order m 4 . For 
example in the second term of eq. (» Eq and (at least for a physical polarization) E\ 
are of order m, while ((2fc 0 ) _1 d/d/co) 2 d^ is of order m 2 due to tree-level improvement. 

In terms of Feynman diagrams, T (1) is represented by 





Correspondingly we write 


pt 1 ) 


p _|_ ps hl _)_ pg h2 _|_ pgh3 _|_ pgh4 _j_ pgh5 

_j_p V 4 I _|_ pV 42 _|_ p\43 _|_ pVKl _|_ j^W2 _|_ pVK3 _j_ pV 3 _|_ pV 5 


where 

r = 2i£ejXYj^](A;,p,g), 

j i 


(143) 

(144) 
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r ghl = o, 
r gh2 : 


k',k" 
r g h3 _ pgh2 




(145) 

(146) 

(147) 


$ Hk+ k ’~ nHp+k"- fc')Gg? 12 2 2) , (ms) 

k Z K z k 

(149) 


pgh4 _ _ 

k',k",k"' 
r g h5 _ pgh4 

T y41 = $ 5 ( k ' + k"-q)Dxy(k , )D pp ,(k")V Xpj (-k',-k",q) 

j k',k” 

x {g 2 _(k',k")V 12yp ,(k t p t k',k") - 2G**'/'V m ,(k,k',p, fc")}, 
r 1 ' 42 = -^E e l $ Kk'+ k"-k)D KX (k')D^(k")V^(-k',-k",k) 

j k',k" 

x {g 2 _(k', k")V j2 \' p '(q,p, k', k") + 2 G k k J pq k V jX > 2p r{q, k',p, k")}, 


p\43 _ pV 42 

T W1 = $ 5(k'+k"-q)D xx ,(k , )D pp ,(k")V Xpj (-k',-k",q) 

j k', k" 


(150) 

(151) 

(152) 


9-{k’X) 
9-(k,p ) 


x{g + (k',k")g + (k,p) + (k'-^k) + (k'^p)} W 12X y(k,p, k', k"), (153) 

Mk'X) 


r^ 2 = bi'E* $ 5(k'+kr-k)DMV)D&(k?)v Xpl {-k',-kr,ky 


j k',k" 


9-(k,p) 


^IV3 


x{g+(k',k")g + (q,p) + (k'<->q) + (k'^p)}W j2X > p >(q,p, k',k"), 

■pW2 


(154) 

(155) 


m = -iE 1 


J 


8{k + k'~ k'")5(p + k"— k')D xx ,{k')D pp ,{k")D TT ,{k'") 


j k',k",k'" 
yk'k'kk'"- 


XGI; P Y y yn{k\ -k m , k)VMk", -k',p)V T , pj (k"\ -k\ q), (156) 

r^ 5 = $ D Xp (k')xk> { [V 12 j Xp {k,p, q, k', -k') + 2 cyclic perms] 

3 k ' 

+ [9-( k ,P + 2 k')V 1X2jp {k , fc',p, q, —k') + 2 cyclic p.]} . (157) 

In the latter equation, the cyclic permutations act on (k, 1), (p, 2) and (q,j)- Furthermore, 
g- is defined in eq. ( |118| ) and 


9 +{k i,fc 2 ) = 2 cos ^(fci,fc 2 )^ , 

fc 4 ,fc 5 ,fc 6 — 9-(k\, kj)g-{k 2 , k 5 )g^(k 3 , fc 6 ) 


G 




^/UM2/43M4 1 1 k 2l k$i k^ — V) 


g-(ki,k 2 ) 

/41/42/43/44 (^1) k 2l k^i k 4 ) ^12/41/43/M 


(158) 

(159) 

(k 2 ,ki,k 3 ,k4). (160) 


The equality of various Feynman diagrams is due to lattice symmetries and properties of 
k,p and g. We checked that in leading order in m the expressions reduce to their infinite 
volume counterparts, cf. the discussion below eq. ( 118 ). Numerically we found that the 
r w contributions are at least of order m 4 . 
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Liischer-Weisz 

square 

L 

N=2 

N=3 

N=2 

N=3 

4 

-0.78341711803619 

-1.47968692446631 

-0.79859872384707 

-1.51787874584542 

6 

-0.99849061664701 

-1.78876943596789 

-1.01681082402891 

-1.81613837864479 

8 

-1.13963010017710 

-1.98424388671692 

-1.15382728105870 

-2.00463419279200 

10 

-1.24091978060661 

-2.12599811096867 

-1.25207004886790 

-2.14258204874697 

12 

-1.31941153232092 

-2.23737552735037 

-1.32855710935494 

-2.25172969936027 

14 

-1.38345219535201 

-2.32928849922146 

-1.39125843561358 

-2.34224049339430 

16 

-1.43757858326173 

-2.40765025540042 

-1.44445943234635 

-2.41966813031578 

18 

-1.48449250619699 

-2.47601715882150 

-1.49071210672691 

-2.48738335434952 

20 

-1.52592324035698 

-2.53669577194007 

-1.53165594879818 

-2.54759007845263 


Table 4: Raw data for X^/m. For each entry the absolute accuracy is estimated to be 
1CT 13 . Note that the inexact formula ( |137| ) was used, so that the data has systematic 
deviations of order m A . 


The analytic expressions for the second and third terms in eq. (|137|) can be found by 
differentiating eqs. (|lll|) -( pT7| ) (the differentiation can be brought over to the integrands 
without problem). We checked the corresponding computer programs against our old 
programs for 7 r^ w (k) by running the latter for near values of k 3 . 

The diagrams proportional to c[ can be calculated analytically. We separate their total 
contribution (which can be read off from eq. ( |136|) ): 


AW 

m 


AW 

m 


+ 4m 2 [9(0^ 


c' 2 ) + 2c' 2 \ + 0(m 4 ). 


(161) 


Like m^j, AW is to be evaluated numerically for a number of lattice sizes. The remarks in 
the previous subsection concerning the integration routine are also valid in the present case. 
Thus the translation of the analytic expressions to computer programs is straightforward. 
The resulting data is listed in table [4j. 

The data is expected to be of the form [|nj (with new coefficients a* and 5;) 


AW 

m 


Oq 


b 0 In m + aim 2 + a 2 m 4 + b 2 m 4 


In m + 03 m 6 + b 3 m 6 In m + 


(162) 


Indeed it can be checked that our data is consistent with the absence of odd powers in m. 
Also we checked the absence of an m 2 lnm term, expected due to tree-level improvement. 
Furthermore, since lim m ^ 0 (A /m) tree i eve i = —8go and A is a renormalizable parameter, b 0 
should equal 

, „ 11N 

b o = 8A) = (163) 

The fit ( |162|) (dropping a 3 or b 3 ) to our data for L — 10 • • • 20 reproduces this value to six 
digits. This a non-trivial check against programming errors (and other errors)^} Moreover, 

16 Useful intermediate checks can be obtained by extracting the coefficients of the m° In m terms for indi¬ 
vidual diagrams, and comparing them to the 1/(4— d) poles of their dimensionally regularized continuum 
counterparts. 
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by fixing b 0 to have the exact value ( |163|) the other coefficients can be fitted to a higher 
accuracy. In this way we obtain 


a 


LW 

0 


a 


LW 

1 


sq 

a 0 


a 


sq 

1 


-0.84832346(3) 

for N = 2 

(164) 

-1.28773532(5) 

for N = 3, 

0.419861(6) for 

N = 2 

(165) 

0.78417(3) for 

N = 3, 

-0.85183887(3) 

for N = 2 

(166) 

-1.29656105(4) 

for N = 3, 

0.324745(5) for 

N = 2 

(167) 

0.59095(2) for 

N = 3, 


where the error-estimate procedure described in the previous subsection was used. 


The LW coefficients agree with those obtained by Liischer and Weisz 


,LW(jy = 


2) = -0.8483231(3), a^ w (A^ = 3) = -1.2877352(1), a™ (N = 2) = 0.41988(3), a\ w (N = 


3) = 0.78412(5). 


"o 
The 


aq coefficient for the square and LW actions differ. This is a 


renormalization effect. Like in subsection |5T| one can see that the a 0 coefficient is related 
to the Lambda parameter: 

( Q o({c*})-ao({ci})) ( 168 ) 


A* 

X 


We thus find 


A 


LW 


A 


sq 


= e°o 

1.0095074(2) 

1.0159636(2) 


for N = 2 
for N = 3, 


(169) 


in complete agreement with eq. (|8T]). Incidentally, for N = 2 we also did a run for the 
Wilson action and found 


a^(N = 2) = -1.37530949(10), (170) 

so that 

4.1308934(14) for N = 2, (171) 

4.0919894(14) for N = 2. (172) 

After renormalization one reads off the improvement condition: 

4 [9(c) - c)) + 2c)] = — 0 , 1 . (173) 

Hence together with the results from the previous subsection we have completely deter¬ 
mined the one-loop Symanzik coefficients for the square (and Liischer-Weisz) action. 

We are convinced of the correctness of our new result ( |167| ) because of the many internal 
checks in our computation. In particular, it is extremely implausible that an expansion 
giving the correct values for the leading coefficients &o and ao (as verified from the Lambda 
ratios) as well as for the subleading coefficient bi, would produce an incorrect value for the 

other subleading coefficient a±. The fact that for the LW action all results by Liischer and 

Weisz have been reproduced, confirms our confidence to a large extent. 


Alw 

A w 

Asq 

Aw 
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7 Summary 

In this paper we computed Lambda ratios, one-loop Symanzik coefficients and the tadpole 
parameter for the Liischer-Weisz and square actions. The results are summarized in table []. 
We remind the reader that coefficient combinations other than C\(g^) = Ci(gg) + 4c 4 (gg) 
and c 2 (go) are unimportant because only Ci(gg) and c 2 (gg) couple to the a 2 corrections of 
any on-shell quantities (see section |j). 



Liischer-Weisz 

square 

Co 

0.135160(13) (N = 2) 
0.23709(6) (N = 3) 

0.113417(11) (N — 2) 
0.19320(4) (N = 3) 

cl 

-0.0139519(8) (N = 2) 
-0.025218(4) (N = 3) 

-0.0112766(7) (N = 2) 
-0.019799(2) (N = 3) 

C 2 

-0.0029431(8) (N = 2) 
-0.004418(4) (N = 3) 

-0.0029005(7) (N = 2) 
-0.004351(2) (N = 3) 

A/ A W il S on 

4.1308935(3) (N = 2) 
5.2921038(3) (N = 3) 

4.0919901(2) (N = 2) 
5.2089503(2) (N = 3) 

/(U) 

0.366262680(2) 

0.3587838551(1) 

/(1.2) 

0.662626785(2) 

0.6542934512(1) 

0 (2, 2) 

1.098143594(2) 

1.0887235337(1) 


Table 5: One-loop improvement coefficients, Lambda parameter ratios and one-loop expec¬ 
tation values of small Wilson loops for the Liischer-Weisz and square Symanzik actions. For 
quantities that were extracted both from the static quark potential and from the twisted 
spectroscopy, we used the most accurate result. In particular we made use of ref. |TB| for 
the combination C\ — c 2 , cf. subsection |5.2| , which we verified to somewhat lower accuracy 
in subsection |6.2| . For the definition of /(L, T) and its relation to the tadpole parameter, 
see eqs. (|34|), (j52|) and 


For the Liischer-Weisz action we found complete agreement with previous calcula¬ 
tions |T5], |1(J. In particular we reproduced the values for the Lambda parameter and 
the one-loop coefficients, as well as for /(1,1) and /(1,2) (which are related to the tad¬ 
pole parameter u 0 )- Especially the agreement with ref. [il5l is non-trivial because we used 
covariant instead of Coulomb gauge fixing. For the square action we computed the com¬ 
bination c\ — c 2 and the Lambda parameter in two independent ways, namely using the 
static quark potential and finite volume spectroscopy. The Lambda parameter was also 
computed in ref. [j3lj from a background field method. The agreement between all results 


leaves us with no doubt that our values for (c\ 


C 2 ) square and A square are correct. We are 


square ? 


also convinced of the correctness of our result for the combination (9(5^ — c 2 ) + 2c' 2 ) 
because it was found from the same numerical data that gave correct results for three 
verifiable coefficients. 

It is interesting that tadpole improvement of the tree-level square action captures 79% 
(SU(3)) to 80% (SU(2)) of the one-loop correction to the appropriate ratio of the coefficients 
Ci(gl) 0, especially because this is similar to the prediction found for the Liischer-Weisz 
action (76% resp. 80%). 
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A Structure of the Feynman rules 


In this appendix we expand the total action, eq. ([H]), in powers of g 0 , largely adopting the 
notation of refs. 0EOI0- We use Fourier space: 


n 2 -i 


A^x) = <?o E / e^+^A^k)^, 
6=1 Jk 


where 


3 

n 


T / a dk„ 


(174) 


(175) 


7* fi—0 27T ) 

It should be noted that the above definitions hold for the infinite-volume case. The modi¬ 
fications for the finite twisted volume are discussed in appendix [Dj. From now on, summa¬ 
tions over SU(iV) and Lorentz indices will be implicit. 


The expansion of the measure and ghost actions can be found in refs. [p~5[,|16|l. For 
completeness we copy the results! 'j: 


N 


1 


^measure = -5 / (2 tt) 8(k x + k 2 ) A^h) A a £(k 2 ) + O(g 0 ), (176) 

CL" Jki,k 2 


S t 


ghost 


[ c ai (k 1 )c a 2 (k 2 ) [(27r) 4 h(/c 1 + k 2 ) 8 aia2 kl + igof ai a 2 a 3 

Jki,k 2 L 

x f (2tt) S(k 3 T k 2 T kAA^(k 3 )ki^ l c 2 ^, T -rrgo 8 fi 3 ii 4 fai a 3 ef a 2 a 4 e 

J k3 1Z 

x a 2 [ (2vr ) 4 8 (h+ k 2 + k 3 + k 4 )A2(k 3 )A^(h)k^ 3 k 2l , 3 + O(g 4 0 ) 

Jk 3 ,k4 


(177) 


Here c“, c a are the Fonrier-transformed ghost fields. Also we adopted the conventional 
notation k 2 = o ^ and 


K = | sin(§*;„), 

= cos(§/c M ). (178) 

The difference between and the action coefficients c, should always be clear from the 
context. 

17 There is a slight overall difference with ref. jl(]] because we explicitly keep a factor (—1/go) in the path 
integral, eq. (p2|), instead of absorbing it into the action. 
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The remaining part of Stotai/<?o is expanded as follows: 


i 00 n n ~ 2 

r(S({c,(9o 2 )}) + S g ,) = E ^S„(Wjo)})- 


9 0 n = 2 

5 n({Q(^)}) is linear in Cj(gg) and therefore we have 

s„( W 9 „ 2 )}) = s n ({c}) + So 2 (S„({c'}) - 5 n ({ 0 })) + O(/ 0 ). 


(179) 


(180) 


In sections |j and H the terms — ^({O})) are treated as insertions. 

Only S 2 , S 3 , S 4 and S$ are needed in our calculations. These can be written in the 
following way: 


Stiici}) 

&({*}) 


= Sr, 


where we defined 


W, 


Iki,k 2 


(2nyS(k 1 + WAZWAZfaXD- 1 )^!), 


Sittci}) = 


fa 1 a 2 a 3 / (27t) 4 <5(/Ci + k 2 + k 3 ) 

Jki,k 2 ,k 3 

xiji(A: 1 )i«(A :2 )i«(A : 3)V ww/1 3(A ! i,A :2 ,A ! 3), 

(27r) 4 5(A;i + k 2 + k 3 + k 4 )A*\(ki)A^(k 2 )A°\(k3)A° 4 (k 4 ) 


Jki,k 2 ,k 3 ,k 4 

^ [faia 2 efa,3a4e(V[iifi 2 H3fi4 (k\, k 2 , /C 3 , /C 4 ) ^/i 2 (Ul/ 43 /i 4 (^ 2 > k\, k 3 , ^, 4 )) 


5 5 ({ci}) = a 


5 , aia2a3a4bf 7 /Ji/i2A t 3m(^'l’ ^2 j k 3 , ^ 4 )] , 

aia2a 3 a 4 a5 / (27r) 4 <f(fci + k 2 + &3 + + ^ 5 ) 

J kiMMMte 


xA^{ki)---A^{h)V, 


M5 v 


^1/22^3^4^5 


(ki,k 2 ,k 3 ,k 4 ,k 5 ), 


V, 


M1/X2M3 




= X>d;> 

= E^ 


AU/72M3 ’ 




^1M2A^3M4 


M1M2M3M4’ 




^1M2M3M4^5 


= v c y(*) 

Z^ L *>if 
i 

= V cT (i) 
Z^ c * >d 


/U/42M3M4’ 


A 1 1/92^13/^4/45 " 


The definitions of the color factors appearing in eqs. ( |183| ), ( |184|) are: 

1 


Sabcd — 

C^abcde 


— Tr (T a T fe T c T d + 23 permutations) , 


rjT^. ^^arpbrpcrpdrpe _ rj^erj^drjncrj^brjna'^ 


(181) 

(182) 

(183) 

(184) 


(185) 


(186) 

(187) 


For our purposes it is not necessary to work out these factors for general values of the 
indices. Note that W /il/i2M3(U4 (A:i, k 2 , k 3 , k 4 ) is completely symmetric under permutations of 
1 • • - 4 and vanishes in the continuum limit. 

The Feynman propagator D^ and its inverse are discussed in appendix while in 
appendix y the vertices are given. The following parametrization, valid as long as only 
planar Wilson loops are included in the action, is convenient: 
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^iiA* 2 w( fc i> fc 2 , A: 3 ) = {h, fa, k 3 )S i 


'M1M2M3 


sfiis^ 1 ’^, &3)<W 2 + 2 cyclic permutations 


Ml M3 

^MlM2M3M4( fc l> fa, fc 3 , fc 4 ) = 4 ? (All, fa, fa, fa)S A^1M2M3M4 

W ik k k t,'iA +3 cyclic perms 


+ 


3^4 (*U, Afe, Afe, ^ 4 )^ 


Ml M2 M3 

(^i’ fai k 3: A: 4 )(5 AllM3 5 M2/Lt4 

+ K ( Ss( fcl ’ fc 2 ’ fc 3> & 4 )<WA 3 M 4 + 1 cyclic perm 

W /_ MiM2M3M4( fc l, fc 3 , fc 4 ) = fjw\ki, k 2 , fa, k 4 )S M1M2M3M4 

■9,S ( fc i) fc 2 , fc 3 , h)S^ lltl2fl3 + 3 cyclic perms 
C^, ^ 2 ’ ^ 3 ’ Wi/Am. + (2 3) + (2 <-► 4) 


14 


M1M2M3M4M5 


(Ah, fa, k 3 , fa, fa) = fj®(k 1 , fa, k 3 , fa, fa)8 t 


■'M1M2M3M4M5 


^MlM5 (^1 ’ ^ 2 ’ ^ 3 ’ ^ 4 ? ^ 5 )<5/11/12M3M4 

h( n% 2 (fa, fa, fa, fa, fa) 8 , 

ht Ufa, fa, fa, fa, fa)^ 


+ 4 cyclic perms 


M 1 M 3 M 5 ^M 2 M 4 + 4 cyclic perms 


MlM2M3 u M4M5 


+ 4 cyclic perms 


(188) 


(189) 


(190) 


(191) 


where 8 ^...^ = (5 Ml/i2 ■ ■ ■ Permutations act simultaneously on Lorentz indices and 

momenta. From considerations given in ref. (T3j it follows that vertices n even, are 

real and invariant under inversion of all momenta. For n odd they are imaginary and odd 
under inversion. Moreover, the components have the following properties with respect to 
permutations: 


fj?\fa, k 2 , k 3 ) 
9^(fa, fa, fa) 
f {4 \fa, k 2 ,fa,fa) 
9$(fa, fa, fa, fa) 
h^ik^h 2 ,k 3 ,k 4 ) 

h'$(fa, fa, fa, fa) 
(fa, fa, fa, fa) 
9 { ™\fa,fa,fa,fa) 
h^\k fa,fa,fa) 

(fa, fa, fa, fa,h) 
9$(fa, fa, fa, fa, fa) 

fa®, (fa, fa, fa, fa, fa) 
h'j$ (A:i, k 2 , k 3 , k 4 , fa) 


-f®(fa,k 2 ,k 1 ) = f®(k 2 ,fa,k 1 ), 

~9$(fa, fa, fa), 

f [ ^( k 4, fa, fa, fa) = f®(fa, fa, fa, fa), 

9$ (fa, fa, fa, fa), 

fa 4 } (fa, fa, fa, fa) = fa® (fa, fa, fa, fa), 
h’$(fa, fa, fa, fa) = h'®(k 3 , fa, fa, fa), 

f^ikpu, k P(i), k P(i), k P{ 1 )) f° r all perms P , 

g^\ k 3 , fa, fa, fa) = g^\fa, fa, fa, fa), 

h^\fa, fa, fa, fa) = fa™\fa, fa, fa, fa) = 

fa$?(fa, fa, fa, fa), 

~f^\fa, fa, fa, fa, fa) = f®(fa, fa, fa, fa, fa), 

-9$(fa, fa, fa, fa, fa), 

-h$(k 5 , fa, fa, fa, fa), 

-h'^(fa, fa, fa, fa, fa). (192) 
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B Propagator 

Following the notation of ref. ||, the inverse propagator as defined in eq. 


reads 


(zrV(*0 = 




J [LV 


q^{k) — 
a 


Ay Ay. 


(193) 


For the action in eq. (f|) the tensor q )iu equals 

Qiivi.k') (1 S^u {k') i 

s fiv{k) = (co + 8 ci + 8 c 2 + I 6 C 4 ) — a 2 (ci — c 2 + 4 c 4 )(p + A: 2 ) 

— a 2 c 2 k 2 + a A cJz 2 Jz 2 ). (194) 

(In view of the discussion in section |3| we do not include C3). Note that, by definition, this 
tensor factorizes for the square action (c 2 = 0, C 0 C 4 = cf). 

The inverse of eq. (|193| ) reads 

1 


D^(k) = 


(k 


2\2 


T,A, f (k)kl 


fifiv [^4/ii/(A^) cr] Ay Ay 


p+v 1 


with A^ satisfying 

A^(k) = 0 (V/i), A^{k) = A Vfl (k) = A^(-k), lim A^ v (k) f = 

a ^° VW) 

In ref. [j5|] the general form of A^ u in terms of q^ can be found. Using the notation 

k {n) = i:K, 

P =0 

fl, v : such that \x 7 ^ v 7 ^ jl 7 ^ z>, 

p — £. 2 1 2 

1 [IV A/ /} A '£o 

Qfii/ — kfj + Ay ; 

we here give the explicit formulas for the Wilson (W), Luscher-Weisz (LW) 
square (sq) cases: 


A^\k) 

A^\k) 


= 1-5, 


fiui 


1-5 


[LV 


(k 2 + h(k (4) + k 2 kfy ( k 2 + ±(k (4) + k 2 kj )) , 


A( lw ) 

A UW) = (p + ^(4)) [fc2 + ^( ( jfe2 ) 2 + j& (4)j 


A^\k) = 


+m (( k 2 ) 3 - PP 4) + 2jfe( 6 >) 
1-5 , 


+ik ^ 2 n ^p> 


(195) 

(196) 

(197) 

(198) 

(199) 

( 200 ) 

and 

( 201 ) 

( 202 ) 

(203) 


J [IV 


{( k 2 + &h (4) ) [P + ^fc (4) 


A( S( i) 

+ T2 _ k(^ 2 ) 2 + I 

+ P< 


P =0 


Q iw + ^k 2 Q 2 


1 n3 


[LV 


Q 


[IV 


T 24 1 


1 Af'p'i 2 _ ip£.( 4 )_Wpi 2 h( 4 ) -i- j-tU 4 )! 2 _i_ 1 

144 m / e tv tv 144 ^ \ 72 IA ) 1 6 

:(p + i*< 4 ) + _|_pfc( 4 ))g Mi/ + ^ (p - fk^)Q% - 
(sM - k 2 + gAP-*) + j 

a^) = (p+ip 4 )) 2 n(i+M 


p=0 


(204) 

(205) 
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In these formulas we used lattice units. The a dependence can be reinstated by substituting 
k —> ak. 


For special momenta and polarizations the propagator greatly simplifies [[L5|,[l6|]- From 
eqs. (|193| ), (|194|) one deduces that vectors e satisfying the following two conditions are 
eigenvectors of (ZW 1 ) ^(k). 


1. The gauge condition E /; k^e^ = 0 is satisfied; 


2. 3p e 1R (or (D) such that V/i e {v\e v ^ 0} k^ e {p, —p}. 

(Solutions to these conditions are for example {e = ( 1 , 0, 0, 0); k — (0, ki, k 2 , k 3 )} and {e = 
(0,1,—1,0); k = (ko,ki,ki,k 3 )}). Moreover, the eigenvalue equals simply J2 P s p.p(k)kp, 
where ji is an arbitrary element of {//| E^ ^ 0}. Therefore, whenever the two conditions are 
satisfied, £ is an eigenvector of the propagator with eigenvalue 


E 


/J.U £ flD (fc) E J, 


= dfi(k), d n (k ) = 


E e 2 

C L 


wave function renormalization Zo(k). They are defined through (cf. subsection [67 

Z 0 (k) 


(206) 


Ep s p,p(k)kp 

In section 0 we need the (tree-level) energy F'o(k) for d\(k), as well as the associated 


di(k) = 


kl + E%( k) 


+ (regular in k 0 ), 


(207) 


for k close to the physical mass shell. An equivalent definition of Z a (k) is 


^o _1 (k) 


1 d 
2 k 0 dk 0 


di\k) 


ko=±iEo(k) 


(208) 


Analytic expressions for ZZ 0 (k) and Z 0 (k) can be easily obtained because d l 1 (k) is quadratic 
in kl for any values of co, Ci, c 2 and C 4 . When cq is fixed by eq. (|23|) the results read 


E 0 (k) 

= 2asinh^/—^ (l - y/l + 4AC/ 5 2 J, 


Z 0 (k) 

- sf'Aw) 

(209) 

where (c 4 = C\ + 4c 4 ) 




A = ci — c 4 fc 2 , 

B = 1 — (ci — c 2 )k\ - 2c 2 k 2 , 

C = (l - c 2 k 2 ) k 2 - (c 4 - c 2 ) (fc 2 k 2 + k (4) ) + c 4 fc 2 k (1) , (210) 


with k ^ 4 - 1 = Ej =1 kj. For c 4 = 0, iZo(k) reduces to the expression found in ref. IlSj , where 
however Zo(k) is incorrect as it has the wrong limit for k —> 0 (which corresponds to the 
continuum limit). The small-fc behavior of Eo(k) and Zo(k) is 


£o(k) 

Z 0 (k) 


1 — l( c i — c 2 + 12 ) ( k" + 


k (T' 

k 2 


+ 0{k 4 


1 — 2 (ci — c 2 + ^)k 2 + (ci — c 2 )kf + 0 (k A ). 


( 211 ) 


Note that unlike in the continuum, Z 0 (k) need not be 1. 
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C Vertex components 


In this appendix we list the vertex components defined in eqs. (|188|) - (|191| ). We restrict 
ourselves to the contributions of Wilson loops with non-zero tree-level coefficients in the 
square action, i.e. the 1 x 1 {i = 0), 1 x 2 {i = 1) and 2 x 2 {i = 4) loops. The results 
for i = 0 and i — 1, as well as for the non-planar cases i = 2 and i = 3 were obtained 


before by Weisz and Wohlert |1(|. We have greatly benefited from comparing to their 
expressions. For our calculation we used Mathematica |26|. We automized the translation 
to only doing by hand line-breaking and some minor changes to improve readability. 

We therefore expect the expressions below to be free of typesetting errors. We are also 
confident that no other mistakes were made because 1) the agreement with Weisz and 
Wohlert’s expressions 18 ; 2) the correctness of our expressions for small lattice spacings a. 
An example of the second point is given at the end of this appendix. 

For shortness we suppress the momentum arguments (k\, ■ ■ ■, k n ) on the left hand sides 
of the formulas below, and set a = 1. The lattice spacing can be reinstated by dimensional 
analysis: the mass dimension of hi equals +1, that of any n- vertex component equals 4 — n. 

Due to the absence of Lorentz symmetry, explicit expressions for the vertices tend to 
be long, especially for i ^ 0 and n > 4. We have put some effort in the search for relatively 
short forms. Nevertheless, we apologize for the lengthy formulas below. 


r(3,i=0) 

J 11 


= 0, 


(3,i=l) _ 


= — i < k{ Ci M (k -2 — k 3 ) + 2 cyclic perms of the momenta >, 


(3,<=4) _ 


(3,<=0) 

y fw 

(3,<=1) 

y fiu 


9 , 


(3,i=4) _ 


fils' 


—i j(2Aq) C\n (fc 2 — k:i) lt + 2 cyclic perms of the momenta [>, 

-*c 3/1 (ki - k 2 ) u , 

-i 14 (cos (k 3/i ) c lfi c 2 fi + cos |(fci - k 2 ) v c 3/1 c 3 „) (k t - k 2 ) l 

-K 2 ^) 3/ , hu(k i - k 2 )Jj, 

-i c 3u {8 cos (k 3fl ) c lfl c 2fl (2(k 1 - k 2 )) u - (2 k) 3fi (2 k) 3u (h - k 2 ) /x \, 


(4,i=0) 


f, 

c(A,i=l) 
J fl 


— y) 4(^1 + ks) — ?(ki + k 2 ) — 2 (ki + ^ 4 ).+ ki p k 2p k 2 p kpp\ , 


( 212 ) 

(213) 

(214) 

(215) 

(216) 

(217) 

(218) 


p=0 

3 


= E 


p=0 


(40 — 2 k 2 lu - 3 k\p -3k 2 4u -3 (h + k 2 )p - 3 (ki + h) u - 2 k 


1/4 


v 2 fi 


M/x 


+ 1(^1 + k 2 ) p 1 — 20 + 2 k\ p + 2 k^ + (^i + k 2 ) p + 2 (Aq + + (ki + k 2 ) 


ip 


+ 3 cyclic perms of the momenta > , 


(219) 


18 Apart from some more or less obvious typographic errors in ref. Pq| , we only disagree on the overall 
sign of the 3-vertex. 
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i?(4,z=4) _ 

J u 


E 

P= 0 


um ip ( 32 - 2 kl - 3 kl - 3 kl -3(h + k 2 ) u — 3(k l + k 4 ) 


U 1M 


v 2fi 


4/x 


+i( 2 (^i + & 2 )) p ( — 16 + 2 kf^ + 2 k 2fl + (&i + k 2 )^ + 2(/ci + /C 4 ) 


+ 3 cyclic perms of the momenta 


(4,i=0) 
y fiu 


(220) 

( 221 ) 


0^=1) 
tf {IIS 


(4,i=4) 

jfiis 


h 


b [C3v (kl - k 2 ) v - Ci„ (k 2 - k 3 ) u - k lu k 2v k 3u J fc 4/t , 

-c 4 „ (cos(£;i - £ 3 )*, ( 2 k) 2u + H 2 ^) 4i ,) & 4 M + cos |(/ci - £ 3 )^ & 2i , 
x (Es(/c 4/i ) (2 c 4/i c 3At k 2ll - c 2fl (k x + /c 3 ) J - 4 c 4/i c 2m c 3m (2/c) 4/i 

12 (2^) 4 ^ k 4l/ ^16 Ci M c 2m c 3 ^ (k 1 k 3 )^ k 2 fi -t~ 2 c 2 ^ ki^ k 3 ^ 

-cos(A; 4m ) c 2/i c 2y (A:i - k 3 )^ (ki - k 3 ) u , (222) 

2 cos(/c 4 - k 3 ) u c 4v ( 2 k) 2u 

x ^cos(A; 4/i ) (2 c 4/i c 3At k 2tl - c 2/i (ki + k 3 )^) - 4 c 4/i c 2m c 3/i ( 2 &) 4/i 

— |c 4y (2A:) 4/i ( 2 fc) 4(/ ^16 c 4/i c 2/i c 3/i — (/c 4 + ^ 3 )^ k 2tl + 2 c 2/i ^ 

-2 cos(k 2u ) cos(fc 4/ J c 2/i c 4i , (fci - fc 3 ) M (2(&i - Afe)),,, (223) 

2 cos §(A ; 2 — ^ 4 )^ cos \(k x — k 3 ) u , (224) 

8 c 4#1 c 3/i cos(/c 2 - fei) M cos |(/ci - £ 3 )^ + 8 c 2u c 4u cos ±(k 2 - /c 4 ) m cos(/ci - /c 3 )^, 

(225) 

32 c 4m c 3At c 2t , c 4v cos(k 2 - k 4 ) fl cos(/ci - A^)*, (226) 

-cos |(fc 3 - h)^ cos |(/ci - k 2 ) v + \k 3ll k 4m k Av k 2vi (227) 

[ - 4 c 4/i c 2/i cos(/c 3 - /c 4 ) M cosi(/c 4 - k 2 ) v - §cos§(A;i + k 2 ) u (ki - k 2 ) fi 

x(2(A* — ^ 4 )) m + c ifi c 2 fi ( 2 k) 3fl ( 2 k) ^ ki u k 2v \ +(ki<-^k 3 , k 2 ^k 4 , fi^u), (228) 
4 (-4 cos(k 3 - k 4 )^ cos(/ci - k 2 ) v + (2 k) 3fi (2k) 4fi (2k) lv (2k) 2t )j c 4m c 2m c 3 „ c 4 „ 

-2 ^c 3l/ c 4i , cos(fci + k 2 ) u (ki - k 2 ) lx ( 2 (k 3 - k 4 ))^ + c 4/i c 2/i cos(/c 3 + /c 4 ) M 

x(/c 3 - k 4 ) u ( 2 (ki - k 2 ))^j—j(ki - k 2 )^( 2 (k 3 + ^ 4)) M (^ 3 - ^ 4 ) 1 / ( 2 (A:i + A: 2 )) v , 

(229) 

f m) = ^ ki p k 2 p k 3 p k 4 p K® (ki, k 2 , k 3 , k 4 ), (230) 


7 (4,i=0) 

U HV 

7 7 (4,i=l) 

u Hi> 


h (4,i=4) 

'V' 


/ (4,i=0) 

/ (4,i=l) 
//.;/ 


,(4,1=4) 

/ij/ 


p=0 


„W») 

y iiu 

h m) 


-hiikivk 2 vk 3 vK^(ki, k 2 , k 3 , k 4 ), 
k 3 iihitkivk 2 uK$(ki, k 2 , k 3 , k A ). 


(231) 

(232) 
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In the last three equations we introduced 


K^= 0) (k 1 ,k 2 ,k 3 ,h) = 2 , 

Kk= l) (ki, k 2 , k 3 , k 4 ) = 32 (ci m c 2m c 3m c 4m + c lv c 2v c 3v c Av ), 
Kuu i.kij k 2l k 3 , /c 4 ) — 512 (c A [ 1 C2fj,c 3 i 1 c A [ i c\i / C2uC 3 uC A u) ■ 


(233) 

(234) 

(235) 


The 5-vertex component f^’ 1 ' is not needed in the calculation of the one-loop Symanzik 
coefficients. The other components are expressed below as polynomials in a small number 
of functions, namely kj\, Cj\ and 

c tj\ = cos h(ki ± kj ) A ; sf jX = i (ki ± kj) x , 

<$jk = cos ( k i ± k j) a; = s ( 2 0* ± k j))\ • ( 236 ) 

For each component only those functions are used that are even or odd with respect to the 
symmetry operation belonging to that component (see eq. (|192|) ). 


vfiv ~ 4ic 5/i |( 3c 23l/ T ^c 23l/ ) s ^ 4l/ 3c 14l/ s 23y |, (237) 

8i |c 5 M C 5 V ((—3c 23j/ + 4c 23 " i ,)s 44l , — 3ffi 4i ,s 23i/ ) 

Tcos(/u 5 ^) ^3 (c 23l/ s 14l/ + c 14i/ s 23 i/) c 23 /i( c i 4 /i + 2— (2c 14/i + 3c^ 4At )c^ 3/i 


(5,i=l) 


+ s 14/i S 23At 


S 14/i S 23/i ) + 2c+ 3v s 14v ((c 14/i + c+ 4fl )(3c 23fl + 5 c 23/j,) 
~ 3s 14/i s 23At + + 2 c 14i/ 3 s 23ai s+ 4m + s 14fl s 23/J ^ 

+ K 2 ^)5 m (6 ( —C 23i/ + C 23u )s l4u (-(c 23/i + c+J S + M + (c 14/i + 


+ ((3^23^ + ^ C 23fi) S Ufi + 3c 14/i S 23At ^ {c23p S 14u + C 14i/ S 23i/) 


+6c 


14i/ 


S 23v (( C 23// + C 23/J S 14m ( C 14/i + C 14/J S 


23/x 


( ( C 23// + C 23/i) S 14/i + ( C 14/i + C 14/J S 23/j) 


23i/ 


(238) 


(5,<=4) 
y [iv 


16 ic 5 u \ cos(fc 5M ) ( 2 c? 4£ ,s^+(- 3s 23/i s+ 4/t + s 14/i ,s+ 3/i ) + 3 (< 


. 2 - 2 - 
23v b Uv 


+ C 


2 — 2 — 
14j/ S 23j/ 


) 


X ( C 23m( C 14^ + 2c 14/i) (2c 14fl + Sc^Jc^ + S 14fl S 23fl $14 fj, s 23 fj) 

+2 c 23u s 14l/ ((c 14/i + c^ 4/i )(3c 23/i + 5 c^ 3/ J — 3s 14/i s 23/i + ^ 

+ ( 2 fc) 5 /i ^6c 14l/ ( —c 23l/ + c 2 +)s 14 „ (— (c 23m + c+>+ 4/i + (c 14/i + 

+ (( 3 c 2 3 /i + 4 c^" 3 At )s 14/i + 3 c 14 /i S 23At ^ ( C 23 i/ C 14 i/ S 14 i/ + Ci+c +,4 


+6 c 44i . 


C 23u S 23u 


« 


C 23 n + C 23^) s 14/i ( C 14 f! + C 14// 


) S 23/i) 


■'23i/ (( C 23^ + c 23/z) S 14ju + ( C 14/i + C 14/J S 23/u) 


’23i/ 


7 (5,i=0) 


-24is 24/1 (c 3l/ cf 5l/ + §£ 3 ^ 15 */) , 


(239) 

(240) 
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ht'? =l) 


7 (5,i=4) 

U HV 


,/(5,i=0) 

u fiu 

i/(5,*=1) 
//// 


i/(5,*=4) 
//// 


- 48?: {(*'24i' + C^J^/z (cos(^)cit + M 2 ^/^ ) + 2 C3/z(c 15/i + 4;Js^ 

X (es^e^ + *fc 3l /sf 5l/ ) + 2c^c 3M sr 5M (c 3l /<4 5l/ + §4/5^) } » ( 241 ) 

— 192i 03 ^( 024 ^ + 0 ^) (s 24/J (c 15/1 + 4 5/i ) + c 24// S 15/z)( COS (^3i/) c l5i/ 4 ' H 2 ^!/ 5 !^) j 


8 * ( c 45// S 45/z + j 

16* {(3c 13/i C2 At + 4 5/i ) C^S^/z + 3C2/z c 1 3v C 2uC. 


2 + _)_ I — ," 2- _|_ r 2 + 
45// w l c 45/z w c 45/t 


) c 45//) 


(242) 

(243) 

’13// 


+ c 45i/( c 45z/+ C 45i/) S 45//+ 2^2 v s 13u ( c 13/z + C U3/z) ( 2 c 2/i s 45/z 2 C 45~/i 4/‘) + C 2/z c 45/z 


XS 


13// 


+ i( c 45// + C 45i/)( 2 ^) 2 ^ S 13!/ S 45// — S 45i/ (§ c 13i/( 2 42i/ + s 4 ^) S^} , (244) 

32* {c 4 + (3 c 13ai C2^ + C 4 5ai ) (c 45l/ + C 45 i/) S 4 5 // 

— 3 c 2 M (c? 3 l/ cos(A; 2 l/ )4k + (c 45/i — c 2 +)c 2 +) (c 45[/ + 4 5 ^)s 13m 

+ f( c 45// + c 45 i/)( 2 ^’) 2 i/ 5 13l/ (c 13/i + C/g^) ^ 2 C 2 ^S 4 ^ — ^ 4 ^^ 2 //) + ^//Qs^S^ 

^45 1 / (3cOs(^)c2//C4 At 5i^5r3 M + (Sc^C^ + C+ ;j ) S^+ S^+ 

+ § c 13//( 2 ^)2i/ ( c 13/z + c 13/z) ( 2c 2//5 4 ^) t — §C45//4 /z) + c 2//C 4 ^ ( S] t 3/i ) | . (245) 


A powerful tool for tracking errors is the check whether the a —> 0 behavior of the 
vertices is consistent with eq. (|23|). As a non-trivial example we perform this analysis for 
the 5-vertex. 

Reinstating the lattice spacing a, the Cj-weighted sums of eqs ( p37|) - (|245|) equal, up to 
0 (a 4 ) corrections, 


9$ = 2ici 2 (c 0 + 20ci + 64 c 4 ) {{k\ — k^) v — 3(fc 2 — k 3 ) u } , 

hjfj = — 12 ia 2 {(co + 20 ci + 64c4)(fc2 — ^ 4 )//+ 4 (ci + 4 c 4 )(A;i —/cs)^} , 

h'^ = 4m 2 {(c 0 + 20ci + 64c 4 )(/c 4 - ^ 5 )//- 6(ci + 4c 4 )(fci -/c 3 ) M } . (246) 


Eq. (P3| ) seems to be violated because h and h! are not proportional to (co + 20ci + 64c 4 ). 
However, we should take into account the remark below eq. (|8|): eq. (p3| ) only holds in 
terms of the held A^, dehned by 


Un(x) = Pexp I / ds A^(x + sfi) 


— 1 + ^ / d m s A^{x + Si/t) • • • A^{x + s m fi). (247) 

m=1 J0<S-t<S2<-< Sm <a 


By expanding eqs. (|247|) , (|T|) with respect to a, the relation between A M and A M is found 
to be 

A^x) = A^x) + J 2 O 2 [\dlA^x') + [AAx'),d^AAx')\] + 0(a 3 ), (248) 


where x' = x + \ajl. Below we neglect the shift over |a/t, because it only contributes total 
derivatives to the Lagrangian density that hence drop out of the action. 
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The commutator term in eq. ( |248| ) is interesting because it mixes n and n + 1 vertices. 
In particular if we expand eq. ( |179| ) with respect to A instead of A, denoting terms by S n , 
we find the relation 


4 JOnl 

5! 55 “ 5! 5s L 


+ 


9o 

4! 


S 4 


'A,,,—. 12 a 2 t erms ^^45 


+ 0 {a 4 


(249) 


It is clear that only the £4 components with a non-vanishing continuum limit, i.e. h^ and 
h contribute to £5 at order a 2 . By comparing the Lorentz-index structures of eqs. (|189|) 
and ( |191|) one easily deduces that only h and h ® are affected by this contribution, 
where of course , h^\ ■ ■ ■ have the same meaning with respect to A as h^\ h^\ 
have with respect to A. After some algebra, including the disentanglement of the SU(iV) 
structure, one finds the explicit result 


gjS = 2m 2 (co + 20ci + 64c 4 ) {(h - k A ) u - 3(k 2 - k 3 ) u } + 0(a 4 ), 

= -4ia 2 (c 0 + 20ci + 64c 4 ) {3(fc 2 - fc 4 ) M + (ki - fc 5 ) M } + 0(a 4 ), 

= 2ia~(co + 20ci + 64c 4 ) {2(fc 4 — fc 5 ) M — (k\ — A; 3 ) M } + 0(a 4 ), (250) 

which is consistent with eq. (]23|) (as it should be). We stress that this is a rather stringent 
test on the structure of and h ( 5, v, and in particular on the signs and numerical 

prefactors of and h'^ 5,l k 


D Feynman rules with a twist 

In a periodic finite volume the Fourier representations given in appendix |A| are still valid, 
as long as one replaces the integral ( |175|) by a sum. In a twisted finite volume the situation 
is different due to color-momentum mixing. The Fourier representation appropriate to the 
‘twisted tube’ (see section jS|j) is given below. 

We use lattice units in this appendix. For shortness we adopt the summation conven¬ 
tion, both for Lorentz indices and for ‘color momentum’ e in the twisted directions 1 and 
2 (e v /m {y = 1, 2) runs over {0,1, • • •, N — 1}, where m = 27t/( NL )). 




N 


9o^ifj-2 


S, 


ghost 


„ S(h + k 2 )A n (h)A K (h){-2^ k y+0(gi), (251) 

ki,k- 

$ c(fci)c(fc 2 


k 1 , ^2 


S(k 1 + k 2 )kf + ig 0 §, S(k 1 + k 2 + k 3 )A l _ l (k 3 )f(ki, k 2 , k 3 
fc 3 

X £ 2 ^ kl,kl )kinC 2 ^ + $ 3(ki + k 2 + k 3 + k4)A /Ji3 (k 3 )A l2 AkA) 


k 3 ,k 4 


&({*}) 

Ssiici}) 


x/(£q, k 3 , e)f(k 2 , h, -e)*-*< c > e >z-h*.*i>fc lM fc 2M + 0(g 0 4 ) 

S(h + k 2 )A^(ki)A^ 2 (k 2 ) (-2 


k 1 , k 2 


5{ki + k 2 + k 3 )A^ 1 (ki)A^ 2 (k 2 )A^ l 3 (k 3 ) 

ki,k2, k 3 

x (— 2 / {ki, k 2 , k 3 )) k 2 , fc 3 ), 


(252) 

(253) 

(254) 


42 













Sid*}) = 


5{ki + k 2 + k 3 + k4)A fll (ki)A fl2 (k2)A fl3 (k 3 )A fli (k 4: ) 


ki,k 2 ,k 3 , &4 


l 


x -2f(k 1 ,k 2 ,e)f(k 3 ,k 4 ,-e)z 2(e,e) 

^ (T)/ 1 (12 /^, 3 /M (^ 1 > ^ 2 ; k 3 , k 4 ) ^ / M2MlA t 3A t 4 (^2 1 ^ 1 ) ^3) ^ 4 ) ) 

- 5(fci, k 2 , k 3 , h)Wfj, lli2 ^ 3fl4 {ki, k 2 , k 3 , h )], 

5 5 ({q}) = $ 5(fci + k 2 + fc 3 + fc 4 + /c 5 )A Mi (/i:i) ■ • -A^ks) 

fell &2, k 3 ,k 4 , fc 5 

Y.C(k\j k 2 , fc 3 , Aq, Aq)^ti/i2A t 3M4A t 5 (^1) ^2; A^ 3 , Aq, Aq). 


(255) 


(256) 


In these equations, and the vertices Id, IF are precisely the same as in infinite volume. 
However, the color factors are now functions of the momenta (in the twisted directions 
only), 


f(h,k 2 ,k 3 ) = —Tv([r kl ,r k 2 ]r k3 ) 


N 


S(ki, k 2 , k 3 , k 4 ) = 


2* (1 - Xfc l+fe+fcs) 

x ^|(fci,fci)+s(fc2.fc2)+s(fci.fc2) S i n ^(jfci, Jfe 2 )^ , 

1 

Tr (r fci r fc2 r fc3 r fc4 + 23 permutations) 


(257) 


24iV 


C(k u k 2 , k 3 , Aq, Aq) = 


s (! - Xfc 1+ fe 2 +fe 3 +fe4) * 2 Ei<*<j<4 (fc - fe d 
x | ;z |(fci+fe,fc3+fc4> cos fc 2 )) cos (^<Aq, Aq>) 

+ (Aq ^ Aq) + (k 2 e-> Aq) | , 

Tr (r fcl r fe2 r fe3 r fc4 r fe5 — r fe5 r fc4 r fc3 r fc3 r fcl ) 

2* (1 - Xk 1+ k 2 +k 3 +k 4 +k 5 ) i<i<i<8 (fci ’ fc ^ 


(258) 


x sin ( 4 E (fc*> fc i) 

JV l<i<j<5 


(259) 


The notation is explained in subsection 0 Note that the factor (— 2z 2( ~ k,k ' > ) appearing 
in S 2 {a} is a color factor, brought about by IV -1 Tr (TfcT-fc) = as opposed to 

Tr (T a T b ) = -\ 8 ab . 


E Positivity of the square action 

For many values of the coefficients Ci(go), the action (|4|) is not positive 0. Such a choice 
of coefficients is unacceptable because the ‘vacuum’ A^ = 0 (i.e. = 1), which has 

zero action, would not be a minimal-action configuration, and hence not a correct held 
configuration to expand about in perturbation theory. It is therefore important to prove 
that the square action is positive. 


43 




The basic ingredient || for the proof is the statement that for any two unitary matrices 
U and V 

ReTr (1 - UV) < 2ReTr (1 -U)+ 2ReTr (1 - V). (260) 

This can be directly applied to eq. (f|) by expressing the 1 x 2, 2 x 2 and corner Wilson 
loops as products of smaller loops (we choose c 3 = 0). For example, one may write 



x 


(261) 


Noting that any unitary matrix U satisfies 

ReTr(l - U) > 0, (262) 

one quickly derives a rough lower bound: 

>5 , ({ c i(fl , o)}) > (co(So) + 8 ( C l(ffo) + 2c 4 (^o)) + T c 2(fi , o)^) “^Wilson, (263) 

where x = min(x, 0). The Wilson action Swiison is positive due to eq. (|262|) . Therefore the 
action (f|) is positive if the prefactor on the right hand side of eq. Q263| ) is. 

At tree level, the Liischer-Weisz || as well as the square action satisfy eq. (|263| ) with 
a positive prefactor. If one includes the one-loop corrections that are computed in this 
paper, table |], the prefactor is still positive for any reasonable value of g%, as long as the 
free parameter c' 4 is chosen not too big. For example, the one-loop square action is positive 
for C 4 = 0; 0 < gl < 10.9, and also for g^ = 1; 144 c 4 < 7.2 (N = 3). 
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